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SECTION  I 
INTRODUCTION 


1.1  PURPOSE 


This  analysis  focuses  on  illuminating  the  logical  and  mathematical 
structure  of  the  location  estimating  algorithms  found  in  the  TRAILBLAZER 
system,  and  identifying  the  assumptions  that  must  hold  for  these  algorithms  to 
give  valid  results.  TRAILBLAZER  is  one  of  several  current  U.  S.  Army  direction¬ 
finding  systems.  These  systems  use  several  lines-of-bearing  to  estimate  the 
location  of  an  enemy  emitter.  Such  a  location  estimate  is  often  called  a 
"fix."  Several  general  methods  for  direction  finding  and  fix  estimations,  some 
with  more  mathematically  rigorous  foundations,  some  frankly  empirical,  are 
discussed  in  Intelligence/Electronic  Warfare  (IEW)  Direction  Finding  and  Fix 
Estimation  Analysis  Report,  Volume  1,  Overview.  The  TRAILRLAZER  algorithms 
analyzed  belong  to  that  most  interesting  hybrid  class  of  empirical  algorithms 
with  a  strong  mathematical  flavor.  Although  the  designer  of  such  an  algorithm 
often  has  a  specific  mathematical  structure  in  mind,  the  empirical  nature  of 
the  algorithm  often  leaves  the  analyst  several  possible  mathematical  interpre¬ 
tations.  This  richness  of  interpretation  increases  the  understanding  of  just 
how  well  the  algorithms  function  in  various  environments  and  how  compatible 
they  are  with  algorithms  found  in  other  systems. 

1 . 2  BACKGROUND 


This  algorithm  analysis  effort  is  being  performed  by  the  Jet  Propul¬ 
sion  Laboratory  for  the  U.  S.  Array  Intelligence  Centex'  and  School  as  a 
research-type  effort  to  increase  the  understanding  of  the  hybrid  mathematical/ 
empirical  algorithms  found  in  intelligence  processing  systems.  Algox-ithm 
results  from  one  system  are  frequently  used  as  input  data  for  another  system. 
Understanding  both  the  assumptions  urder  which  the  algorithms  work,  and  the 
aao  uu«p  ticns  thcxr  results  aohrefy,  ic  crucial  to  unde  r  e  t  and  mg  the  overall 
system.  This  view  of  a  metasystem  of  intelligence  processing  systems  (see 
Figure  1-1)  is  central  to  this  algorithm  analysis  effort. 
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For  purposes  of  these  studies,  "algorithm"  means  a  set  of  rules  for 
carrying  out  a  single  conceptual  operation  on  a  set  of  data.  There  are  many 
types  of  algorithms  necessary  to  the  operation  of  the  metasystem  shown  in 
Figure  1-1.  Analyses  reported  on  so  far,  listed  in  Appendix  E,  have  focused 
on  four  of  these:  geographical  transformation  algorithms,  self  and  cross¬ 
correlation  algorithms,  and  aggregation  algorithms.  Geographical  transforma¬ 
tion  algorithms  translate  locations  from  one  grid  reference  system  to  another. 
These  algorithms  appear  in  almost  all  systems,  often  as  incoming  data  or  report 
preparation  functions.  Self-correlation  algorithms  test  if  the  entity  referred 
to  in  a  new  report  has  already  been  recorded  in  the  database  that  reflects  the 
estimated  enemy  situation.  Cross-correlation  algorithms  test  if  a  sighted 
piece  of  equipment  belongs  to  an  already  identified  unit,  or  a  lower  echelon 
unit  to  a  higher  echelon  one.  Aggregation  algorithms  try  to  identify  an 
artillery  battery  in  a  cluster  of  equipment,  a  division  in  a  group  of 
regiments,  or  like  groupings.  Several  statistical  issues  arising  particularly 
in  the  correlation  algorithms,  are  analyzed  in  a  companion  set  of  technical 
memo  randa . 

Looking  once  more  at  Figure  1-1,  note  that  the  same  intelligence 
function,  hence  algorithms  performing  that  function,  is  often  embedded  in 
several  intelligence  processing  systems.  Some  generic  algorithms,  such  as  the 
geographical  transformation  algorithms  mentioned  above,  appear  in  almost  all 
systems.  Comparing  these  algorithms  that  perform  the  same  function  in 
different  systems  increases  the  understanding  not  only  of  what  these  algorithms 
actually  do  and  how  well  they  perform,  but  also  increases  the  understanding  of 
how  a  "good"  algorithm  would  work  and  what  it  would  look  like.  Such  compari¬ 
sons  should  lead  to  developing  criteria  for  selecting  algorithms  for  embedding 
in  new  or  upgraded  systems,  and  finally  in  the  creation  of  a  library  of  "good" 
algorithms  from  which  the  choice  can  be  made.  The  development  of  these 
criteria  and  building  such  a  library  are  two  major  goals  of  this  algorithm 
analysis  effort  to  which  each  analysis  of  an  algorithm  in  an  existing  system 
contributes . 


SECTION  2 

ASSUMPTIONS,  RESTRICTIONS,  SCOPE 


2.1  BRIEF  DESCRIPTION  OF  RADIO  DIRECTION-FINDING  AND  POSITION  FIXING 

The  purpose  of  radio  direction-finding  is  to  estimate  or  fix  the 
position  of  selected  emitters.  Usually,  the  position  estimate  is  accompanied 
by  a  confidence  region  reflecting  measurement  errors,  propagation  errors,  and 
modeling  errors. 

Radio  direction-finding  (DF)  requires  that  an  emitter  be  viewed  from 
at  least  two  DF  stations  spaced  far  enough  apart  that  their  look  angles  inter¬ 
sect  as  close  to  90  degrees  as  possible.  However,  90°  is  usually  impossible 
under  battlefield  conditions.  Figure  2-1  illustrates  a  simple  situation  of  two 
DF  stations. 

The  fix  estimate  is  at  the  point  of  intersection  of  the  two  lines- 
of-bearing  (LOBs)  (Figure  2-1).  Since  there  is  only  one  point  of  intersection, 
we  have  insufficient  information  to  estimate  the  fix  uncertainty  due  to 
measurement,  propagation,  and  modeling  errors. 

In  a  multiple  DF  station  configuration,  there  are  many  intersections 
(Figure  2-2).  A  more  accurate  fix  estimate  may  be  obtained  by  evaluating  the 
clustering  of  these  intersections.  Since  each  intersection  is  a  simple  fix 
estimate,  the  uncertainty  can  then  be  expressed  as  a  confidence  region  sur¬ 
rounding  this  fix  estimate.  This  uncertainty  reflects: 

(1)  Random  measurement  errors  in  measuring  the  lines-of-bearing. 

(2)  Errors  because  of  different  radio  propagation  effects  along 
the  lines-of-bearing. 

(3)  Errors  because  of  spherical  or  flat-Earth  assumptions. 

(4)  Phantom  or  ghost  intersections  because  of  the  presence  of 
multiple  emitters  or  hidden  emitter  reflectors. 
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Figure  2-1.  Relative  Location  of  Two  Direction-Finding  Stations 


DF  STATION  1 


DF  STATION  2 


Figure  2-2.  Relative  Location  of  Several  Direction-Finding  Stations 
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ASSUMPTIONS 


Some  standard  assumptions  are  made  in  the  following  analyses: 

(1)  The  lines-of -bearing  are  straight. 

(2)  The  errors  in  the  separate  lines -of -bearing  are  independent. 

(3)  The  errors  in  the  lines-of-bearing  are  Normally  (Gaussian) 
distributed  with  zero  mean  and  fixed  estimable  variance. 

(4)  The  emitter  location  estimate  error  is  distributed  as  a 
bivariate  Normal  distribution. 

(5)  The  sensor  positions  are  known  exactly. 

(6)  The  transmitter  location  is  fixed  during  the  period  of  DF 
fixing. 

(7)  The  sensors  are  properly  sited,  calibrated,  and  operated. 

Assumption  #1  is  reasonable  for  the  systems  considered  in  this 
report  when  the  sensors  are  properly  sited.  However,  this  assumption  is  weak 
at  frequencies  below  approximately  30  MHz  because  of  the  effects  of  atmospheric 
tilt . 


Assumption  #2  is  reasonable  based  on  the  systematic  errors  being 
accounted  for  in  calibrations.  This  assumption  is  weak  at  frequencies  below 
approximately  30  MHz  when  some  stations  are  close  enough  to  each  other  to  be 
subjected  to  the  same  propagation  effects. 

Assumption  #3  is  usual  when  considering  measurements  which  are  sub¬ 
ject  to  random  measurement  error.  There  are  biases  in  the  measurements  from 
navigation  errors,  errors  in  the  calibration  tables,  interference,  depression 
angle  effects,  etc;  these  biases  may  be  removed.  In  the  absence  of  specific 
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knowledge  about  these  errors  tne  normal  assumption  is  reasonable.  Distorting 
effects  such  as  plinthing  to  account  for  wild  bearings,  skewedness  because  of 
low  receiver  signal-to-noi.  e  ratios,  and  distortions  resulting  from  the  sensors 
not  uniformly  surrounding  the  emitter  can  weaken  or  invalidate  this  assumption. 

Assumption  #4  is  necessary  to  allow  confidence  levels  about  the 
estimated  emitter  position  to  be  computed.  The  qualifications  on  assumption  #3 
also  apply  to  #4. 

Assumption  #5  is  reasonable  based  on  the  fact  that  any  such  position 
errors  can  be  added  to  the  emitter  estimate  uncertainty,  if  they  are  signif¬ 
icant  . 


Assumption  it 6  is  necessary  to  the  analyses  of  the  systems  considered 
in  this  report,  and  it  is  reasonable  over  the  period  required  to  obtain  a 
single  fix. 
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Assumption  it7  is  reasonable  in  the  absence  of  contradictory  infor¬ 
mation. 

2.3  RESTRICTIONS 

In  addition  to  the  assumptions  discussed  in  section  2.2,  this  report 
does  not  consider  the  following  effects; 

(1)  Geographic  transformation,  map  projection  effects,  and  grid 
reference  system  conversions  (see  UAA002  Analysis  of  Geo¬ 
graphic  Transformation  Algorithms  July  9,  1982  of  this  series 
of  algorithm  analysis  reports). 


Propagation  effects. 

Centroid  effects  and  susceptibility  to  deception  (meaconing, 
gated  signal  parameter  techniques,  etc.). 


(2) 

(3) 
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(4)  Special  problems  associated  with  low-probabi] ity-of-interccpt 
emitters  (low  SNR,  spread- spectrum,  time-frequency  diversity, 
frequency  agility,  etc). 

(5)  Numerical  computation  and  normal  truncation  effects. 

(6)  Combination  of  lines-of-bearing,  or  emitter  location  estimates 
and  their  confidence  ellipses  from  different  systems  (these 
problems  will  be  the  subject  of  a  future  report  in  this  series 
of  algorithm  analysis  reports). 

(7)  Elimination  of  wild  bearings  and  ghost  intersections  using 
hardware/software  processing  of  target  message  internals. 

2.4  SCOPE 

This  report  covers  the  TRAILBLAZER  system  as  documented  in  ROlM 
1602  Extended  Assembly  Language  listings,  marked  DSOrTBSYS.SV,  generated  on 
2/25/82. 
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TRAILBLAZER  DF  FIX  ESTIMATION 
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3.1  A  GENERAL  DESCRIPTION  OF  TRAILBLAZER  DF  FIX  ESTIMATION 

This  analysis  of  the  TRAILBLAZER.  DF  Fixing  System  (a  communication 
intelligence  collection  system  (COMINT))  is  based  on  the  TRAILBLAZER  AN/TSQ- 
114,  Operator's  Manual  TK  32-5811-022-10-1,  and  assembly  language  listings 
dated  2/25/82.  There  are  some  questions  as  to  whether  the  manual  and  the 
listing  correspond  to  the  same  version  of  TRAILBLAZER,  of  which  there  are 
several . 

TRAILBLAZER  is  a  ground-based,  computer-assisted  COMINI’  DF  Fixing 
System  consisting  of  five  sensors:  two  master  control  stations  (MCC)  and  three 
remote  slave  stations  (RSS).  TRAILBLAZER  can  obtain  relatively  accurate  fixes 
in  the  "normal  fix  mode"  with  as  few  as  three  operational  sensors.  Less  reli¬ 
able  fixes  (cuts)  can  be  obtained  in  the  "degraded  fix  mode"  using  multiple 
lines-of-bearing  from  each  of  only  two  sensors. 

Figure  3-1  depicts  the  most  desirable  siting  of  the  five  stations 
of  a  TRAILBLAZER  system.  This  layout  allows  for  a  maximum  DF  base  line  con¬ 
sistent  with  maintaining  the  required  data-links  between  the  sensors.  Over 
flat  terrain  the  penetration  of  the  system  is  about  15  to  20  km. 

Figure  3-2  indicates  extended  penetration  ranges  possible  when  the 
system  is  operated  from  elevated  vantage  points. 
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TRAILBLAZER ' s  five  sensors  can  obtain  up  to  five  lines-of-baaring 
(LOBs)  simultaneously  on  a  desired  emitter(s)  and  place  this  set  of  LOBs  in  one 
of  up  to  five  available  bins  (arrays).  Each  bin  may  contain  up  to  five  sets  of 
LOBs  on  same  or  different  emitters.  Wild  LOBs  may  be  edited  (rejected)  by  one 
of  the  syctcs  opcjrotorc  boood  on  net vm  1  content. 

The  system  operates  in  a  multifix  (automatic)  mode  and  in  a  single 
fix  (manual)  mode.  Since  the  single  fix  mode  amounts  to  the  first  pass  in  the 
multifix  mode  it  will  not  be  discussed  separately. 
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Figure  3-L.  TRAILBLAZER  Deployment 

From:  U.S.  Army  Field  Manual  (FM  30-476), 

Radio  Direction-Finding,  8  April  1977 
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Figure  3-2.  Normal  Radio  Range  vs  Antenna  Height 


From:  U.S.  Army  Field  Manual  (FM  30-476), 

Radio  Direction-Finding,  8  April  1977 
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3.2  SIMPLIFIED  DESCRIPTION  OF  TRAILBLAZER  DF  FIXING  PROCESS 

The  following  description  is  based  on  obtaining  a  DF  fix  using  one 
set  of  five  LOBs  in  one  bin  (array).  Since  there  are  five  LOBs ,  there  can  be 
up  to  10  possible  intersections  of  two  LOBs. 


rn  _  n] 
r  ”  r!  (n-r)! 


is  the  number  of  combinations  of  n  things  taking  r  at  a  time. 

Each  of  these  intersections  represents  an  initial  fix  candidate. 

The  "best"  intersection  must  be  selected,  and  any  obviously  ".wild"  (extraneous) 
LOBs  must  be  discarded  (edited).  Figure  3-3  represents  a  set  of  five  LOBs  with 
the  intersections  numbered  for  reference. 


The  first  step  is  to  edit  any  wild  bearings.  Since  LOB  from  station 
five  docs  not  form  any  intersect  ions  near  the  cluster  of  the  ouhei  iuL£i’Sec~ 
tions  (see  Figure  3-3),  it  will  be  edited  from  the  set  of  LOBs  as  a  wild 
bearing  (Figure  3-4).  "Ghost"  or  "phantom"  intersections  from  the  geometry  of 
the  LOBs  should  also  be  edited  (Figure  3-5).  These  are  inadvertent  crossings 
of  LOBs  and  not  relevent  to  the  fix  estimation  process. 

Next,  each  of  the  remaining  intersections  is  evaluated  to  choose  the 
"best"  one.  The  "best"  intersection  is  determined  by  considering  which  inter¬ 
section  is  best  supported  by  the  other  LOBs.  For  a  LOB  to  support  a  fix  esti¬ 
mate  the  "exactness  of  the  LOB"  from  the  station  to  the  estimated  fix  is  deter¬ 
mined.  If  the  angular  difference  between  the  two  LOBs  is  excessive  (greater 
than  3  standard  deviation  (sigma)  units  in  statistical  terms*)  the  station's 
LOB  is  considered  to  be  a  wild  bearing  and  is  discarded.  Otherwise,  it  is  con¬ 
sidered  a  supporting  LOB.  The  number  of  supporting  LOBs  is  noted  for  each 
intersection.  The  intersection  with  the  greatest  number  of  supporting  LOBs  is 
selected  as  the  initial  fix  estimate.  Figures  3-6  and  3-7  show  how  the 
$ijnnnrtinff  LOBs  ^re  determined . 


*+  one  sigma  about  the  average  value  of  the  LOBs  will  usually  contain  about 
68%  of  all  the  LOBs.  +  3  sigma  corresponds  to  about  99.7%. 


it  i 


Figure  3-3.  Set  of  Five  LOBs  Showing  Intersections 


Figure  3-4,  Set  of  LOBs  Showing  Intersections  After  LOB5  has  been  Removed 


INTERSECTION  1-2 
EVALUATED 
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The  initial  fix  estimate  is  then  optimized  to  obtain  the  "best"  fix 
estimate.  This  optimization  process  uses  a  potential  function  (derailed  belov) 
as  the  objective  function  and  is  based  on  displacing  the  initial  fix  estimate 
systematically  in  four  directions  (north,  west,  south,  east),  and  calculating 
the  support  for  each  displaced  (trial)  estimate  and  keeping  the  "best"  one. 

This  process  is  continued  with  steps  decreasing  in  size  until  a  "best"  fix 
estimate  is  located. 

The  support  of  the  trial  "best"  fix  is  calculated  as  the  sum  of  the 
potential  function  weighted  (which  is  a  semi-normalized  Gauss ian-we ighted ) 
miss-angles  between  the  actual  LOBs  and  the  computed  LOBs  ro  the  trial  loca¬ 
tion.  The  effect  of  the  Gaussian  weighting  is  to  give  more  consideration  to 
the  LOBs  associated  with  the  smaller  miss-angles. 

Having  found  an  optimized  or  "best"  fix  estimate,  a  "confidence" 
region  is  calculated.  A  "confidence"  region  is  a  region  that  is  likely  to 
contain  the  true  emitter  location  for  some  percentage  of  all  fixes  on  a  given 
emitter  (50%  in  the  case  of  TRAILBLAZER ) . 

The  following  description  of  TRAILBLAZER  DF  Fixing  has  been 
freely-adapted  from  the  software  comments.  Differences  between  the  actual 
code  and  the  description  of  the  algorithm  will  follow,  along  with  comments  m 
the  methods  used. 

The  TRAILBLAZER  DF  Fixing  Algorithms  are  based  mainly  on  heuristic 
and  empirical  reasoning,  rather  than  purely  mathematical/statistical  tech¬ 
niques.  The  four  main  steps  in  the  DF  Fixing  process  are  quite  intermingled 
and  relate  to  the  following  discussion  as  follows: 

(1)  Obtain  initial  fix  estimates. 

(2)  Reject  wild  1 ines-of-bear ing  by  manual  screen  editing  and 
automatic  rejection. 

(3)  Refine  (optimize)  the  initial  fix  estimate. 
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(4)  Establish  a  confidence  region  (elliptical  error  probable 
(EEP))  around  the  fix  estimate. 

3.3  DETAILED  DESCRIPTION  OF  TRAILBLAZER  FIXING 

The  TRAILBLAZER  fix  algorithm  has  been  tailored  specifically  for 
operation  with  the  LOB  data  produced  by  a  ground-based  DF  network,  consisting 
of  a  predetermined  number  of  DF  stations,  whose  locations  remain  invariant 
during  the  data  collection  process.  The  goal  of  this  algorithm  is  its  attempt 
to  resolve  multiple  targets  reliably,  while  at  the  same  time  avoiding  ghosts, 
i.e.,  false  targets  arising  from  coincidental  intersections  of  unrelated  LOBs. 


The  TRAILBLAZER  fix  algorithm  execution  consists  of  three  distinct 
processes  or  phases: 


I R  • 
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(1)  Fix  estimation  (ESTMP  Procedure)  and 

wild  bearing  rejection  (ESTMP  and  FINAL  Procedures). 

(2)  Fix  optimization  (FPEAK  Procedure). 

(3)  Computation  of  an  error  ellipse  surrounding  the 
established  fix  point  (FINAL  Procedure). 


r.j 


Figure  3-8  outlines  the  flow  of  the  TRAILBLAZER  DF  fixing  algorithm. 

3.4  FIX  ESTIMATION 

The  estimation  process  is  the  key  to  the  fix  algorithm.  For  an 
unde*  nding  of  this  process,  the  TRAILBLAZER  LOB  database  structure  must  be 
explain^  ...  The  LOB  uata  are  stored  in  sets  of  up  to  five  LOBs,  i.e.,  one  LOB 
from  each  of  five  possible  DF  stations.  A  set  of  five  LOBs  results  from  a 
system  response  to  a  DF  command  (or  from  a  single  manual  LOB  entry  sequence  via 
the  "demo"  command).  The  assumption  is  made  that  to  the  best  of  the  operator's 
judgment,  the  LOBs  within  a  set  are  associated  with  a  single  emitter.  Given 
that  the  LOBs  within  a  set  are  collected  simultaneously,  this  is  a  fair  assump¬ 
tion  . 
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The  objective  of  the  estimation  process  is  to  examine  all  possible 
intersections  of  pairs  of  LOBs  in  the  same  set  and  all  other  sets,  and  to 
select  that  intersection  which  has  the  largest  number  of  other  LOBs  in  the 
database  that  miss  this  intersection  by  less  than  some  specified  miss  angle. 
The  estimation  process  selects  that  intersection  which  lies  within  the  larges'- 
or  strongest  cluster  of  intersecting  LOBs.  It  also  makes  a  gross  check  to 
prevent  duplication  of  previous  fixes  in  multifix  processing  (as  discussed  in 
intersection  criterion  (2)  discussed  below).  This  process  is  computationally 
efficient  because: 

(1)  The  number  of  DF  stations  is  small.  Therefore,  the 
number  of  intersections  to  be  examined  per  set  is 
reasonable . 

(2)  The  locations  of  the  DF  stations  are  fixed.  Conse¬ 
quently,  for  miss  angle  computations  at  each  inter¬ 
section,  there  are  at  most  only  three  other  exact  LOBs 
that  need  to  be  computed,  since  two  are  already  used 
for  the  intersection. 

Each  possible  intersection  of  two  LOBs  is  calculated  directly  from 
the  geometry  indicated  in  Figure  3-9. 

The  calculate,  intersection  is  validity-checked  by: 

(1)  Verifying  a  true  intersection,  i.e.,  and  both 
positive . 

(2)  Verifying  that  the  angular  difference  between  the  LOBs 
is  sufficient  LOB^  -  LOB^  >  0.6  degree. 

(3)  Verifying  that  the  minimum  and  maximum  range  limitations 
are  not  exceeded  for  either  station,  0.5  <  D,.  <  100  km. 

(4)  Verifying  that  the  intersection  does  not  fall  within  the 
error  ellipse  of  any  previously  computed  fix  (in  the 
multifix  mode  only). 
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Supporting  LOBs  are  calculated  from  all  the  available  (unused  in 
any  previous  fix  optimization)  LOBs  over  all  the  bins.  These  supporting  LOBs 
are  calculated  as  indicated  in  Figure  3-10  and  must  fall  within  3  sigma  of  the 
actual  intersection  to  be  considered  as  supporting  it.  The  actual  value  of 
sigma  is  confused  in  the  available  source  code  listings.  It  is  stated  to  be 
2  degrees  in  the  source  code  comments,  but  assigned  a  value  of  8  degrees  in 
parts  of  the  source  code. 

In  order  for  an  intersection  to  be  a  candidate  for  valid  fix  esti¬ 
mate,  an  intersection  must  satisfy  the  following  criteria: 

(1)  The  intersection  is  real,  i.e.,  the  absolute  value  of 
the  difference  oetween  the  LOBs  is  at  least  0.6  degree 
and  the  directed  LOB  vectors  must  intersect.  The  source 
code  does  not  verify  that  only  forward-looking  LOBs  are 
considered  for  intersections.  It  is  possible,  in 
running  the  program,  to  allow  reciprocal  bearings  and 
create  "phantom"  intersections.  This,  however,  is 
unlikely  because  of  the  normal  deployment  geometry. 

Also,  the  intersection  must  be  in  the  range  of  0.5  to 
100  km  of  the  reporting  stations.  The  source  code 
implementation  of  intersection  out  of  range  fails  and 
would  loop  infinitely  because  of  an  initialization 
problem. 

(2)  If  this  is  not  the  first  pass  for  fix  processing,  i.e., 
a  multifix  (as  opposed  to  single  fix)  situation  looking 
for  multiple  emitters,  the  intersection  must  fall  out¬ 
side  the  error  ellipse  of  the  immediately  previous  suc¬ 
cessful  fix  estimate  (this  is  because  all  previous  fix 
error  ellipses  are  checked  in  the  optimization  process). 
This  provides  the  capacity  to  resolve  multiple  targets. 

(3. a)  In  the  normal  fix  mode  (as  opposed  to  degraded)  the 
intersection  must  be  supported  by  LOBs  from  at  least 
three  stations.  The  supporting  LOBs  must  be  within  a 
+3  sigma  miss  angle  of  the  three  exact  LOBs  from  their 
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Figure  3-10.  Supporting  LOB  Geometry 

stations  to  the  intersection.  Also  the  supporting  LOBs 
belong  to  a  set  of  LOBs  not  used  in  any  previous  fix 
optimization  computation,  and  in  which  the  majority  of 
LOBs  conform  to  this  miss  angle  requirement. 

(3.b)  In  the  degraded  fix  mode,  the  intersection  must  be 

supported  by  more  than  one  set  of  intersecting  LOBs  from 
two  stations  whose  LOBs  were  used  in  the  computation  of 
the  intersection.  The  supporting  LOBs  must  be  within 
+3  sigma  of  the  exact  LOBs  from  their  stations  to  the 
intersection.  Furthermore,  these  LOBs  must  belong  to  a 
set  of  LOBs  not  used  in  any  previous  fix  computations  in 
which  the  majority  of  LOBs  conform  to  this  miss  angle 
requirement.  The  degraded  mode  is  used  only  if  the 
estimation  process  for  the  normal  fix  fails  to  yield  a 
valid  estimate  which  would  satisfy  the  three  station 
criteria. 
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(5)  Of  all  the  possible  intersections,  this  intersection  is 
supported  by  the  largest  number  of  LOBs  that  meet  the 
miss  angle  and  set  requirements  mentioned  in  the  above 
criteria. 

The  estimation  process  yields  a  "non-ghost"  intersection  which  is 
the  best  fix  estimate  (lies  within  the  strongest  clustering  of  LOBs).  Also  in 
the  multiple  fix  case,  the  best  fixed  estimate  lies  outside  the  error  ellipse 
of  the  previous  fix.  In  addition,  by  discarding  invalid  intersections,  the 
amount  of  computations  in  subsequent  estimations  (in  the  multitude  fix  case)  is 
considerably  reduced.  Thus,  for  the  multiple  fix  case,  the  number  of  computa¬ 
tions  is  reasonably  bounded. 

In  the  case  of  multiple  target  examination,  once  the  fix  mode  has 
been  degraded,  it  remains  degraded  for  all  subsequent  passes.  Also,  failure 
to  obtain  a  fix  estimate  in  the  degraded  mode  precludes  further  passes. 

3.5  FIX  OPTIMIZATION 

The  fix  optimization  process  seeks  to  improve  the  fix  estimate  by 
finding  that  location  which  locally  maximizes  a  multipeaked  objective  function. 
By  the  nature  of  the  estimation  process,  the  initial  estimate  should  be  fairly 
close  to  the  optimum  fix  location.  This  process  also  performs  a  final  check  to 
eliminate  the  duplication  of  previous  fixes  in  the  multifix  situation. 

At  the  outset  of  the  optimization  process,  an  estimate  is  available 
along  with  the  exact  LOBs  associated  with  it  from  each  system  station.  Also 
available  are  the  LOBs  supporting  this  estimate  as  a  valid  fix  estimate. 

Before  beginning  optimization,  that  portion  of  the  database  which  supports  a 
particular  fix  estimate  is  modified  to  include  complete  sets  when  a  majority 
of  the  LOBs  supported  the  estimate.  This  increases  the  potential  number  of 
peaks  and  ridges  in  the  objective  function. 

The  optimization  is  then  performed  by  systematically  searching  for 
a  local  peak  in  the  total  objective  function  which  is  a  potential  function  of 
the  form 
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where  the  outer  sum  is  over  the  stations,  the  inner  sum  over  the  LOBs  from 
each  station  and  a„  is  the  miss  angle  between  the  actual  and  computed  LOB 
(see  Figure  5-8).  This  objective  function  is  applied  locally  in  each  case  by: 
(1)  using  the  selected  portion  of  the  database,  (2)  using  the  fix  estimate  as 
an  initial  reference  location,  (3)  using  its  associated  set  of  exact  LOBs,  and 
(4)  using  the  computed  potential  function  for  the  point. 


The  peak  searching  scheme  is  a  fairly  conventional  pattern  search 
optimization  method  (Jacoby,  1972;  Gill,  1980)  but  without  a  pattern  step 
directed  along  the  steepest  gradient.  The  reference  location  is  displaced  by 
some  step  size  (initially  16  screen  raster  (resolution)  points)  along  the  axes 
in  the  following  four  directions:  +Y,+X,-Y,-X.  Only  one  direction  is  con¬ 
sidered  at  a  time,  generating  a  trial  location.  If  one  of  these  trial  loca¬ 
tions  yields  a  higher  value  of  the  potential  function,  the  trial  location 
becomes  the  reference,  an  associated  set  of  exact  LOBs  is  determined,  and  a 
new  trial  location  in  the  same  direction  is  attempted.  This  peak  search 
scheme  continues  until  no  further  improvement  in  the  same  direction  can  be  made 
and  the  process  has  been  repeated  in  all  directions.  Figure  3-11,  Parts  I 
through  IV  illustrate  the  fix  estimate  optimization  process.  At  this  point, 
the  step  size  is  halved  and  the  four  directions  are  tried  again.  The  procedure 
stops  either  when  the  step  size  becomes  too  small  (currently  less  than  one 


raster  point  in  screen  geometry),  or  af tei  a  maximum  uimibei  of  successful  im¬ 


provements  have  been  made  (currently  16).  Under  these  restrictions  the  process 


is  nondivergent ,  that  is,  it  stops. 


Although  this  algorithm  always  stops,  and  gives  a  value,  it  may  not 
converge  in  the  sense  that  the  final  location  estimate  is  substantially  closer 
to  the  location  of  the  nearest  local  peak  than  was  the  starting  value.  One 
major  cause  is  that  the  objective  function  itself  may  not  be  sharply  peaked, 
and  may  even  have  ridges.  Thus,  since  the  optimization  algorithm  does  not  have 
a  pattern  step,  and  it  does  not  rotate  axes  to  take  advantage  of  the  gradient, 
it  may  climb  very  slowly.  Three  other  factors  compound  this  behavior. 


(1)  Angles  are  rounded  to  quarter  degrees. 
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Figure  3-11.  Fix  Estimate  Optimization,  Part  I 
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Figure  3-11.  Fix  Estimate  Optimization,  Part  II 
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Figure  3-L1.  Fix  Estimate  Optimization,  Part  III 
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Figure  3-11.  Fix  Estimate  Optimization,  Part  IV 


(2)  Only  the  first  four  terras  of  the  Taylor  expansion  for  the 
exponential  function  are  used  in  calculating  the  objective 
function. 

(3)  The  algorithm  terminates  after,  at  roost,  16  steps. 

Therefore,  as  often  occurs  in  purely  empirical  methods,  the  location  estimate 
may  be  fairly  inaccurate  and  thus  not  well-defined. 

In  the  multifix  situation  (multiple  passes),  the  final  location  (a 
result  of  peak  searching),  is  checked  to  determine  whether  the  emitter  lies 
outside  the  error  ellipses  of  all  previous  fixes.  If  not,  this  optimized  fix 
is  unacceptable  and  a  new  fix  estimate  must  be  obtained.  Mote  that  if  this 
occurs,  the  selected  portion  of  the  local  database  that  gave  rise  to  this 
unacceptable  optimized  estimate  is  removed  from  further  consideration. 

If  the  final  location  is  acceptable,  what  remains  is  the  computation 
of  che  parameters  in  an  error  ellipse  surrounding  the  optimized  fix  estimate. 

The  data  available  at  this  point  is:  (1)  a  selected  subset  of  the  LOB  database, 
(2)  a  final  optimized  fix  location,  (3)  a  set  of  exact  LOBs  to  that  location 
from  all  system  stations,  and  (^)  and  the  value  of  the  potential  function  for 
that  location.  This  procedure  is  satisfactory  for  ellipses  with  small  eccen¬ 
tricity.  However,  it  degrades  with  the  higher  eccentricity  ellipses  arising 
from  emitter  ranges  that  are  large  with  respect  to  the  sensors'  baseline. 

Also,  just  as  it  is  unclear  that  the  best  fix  estimate  is  a  measure  of  central 
tendency  for  a  known  bivariate  distribution,  it  is  equally  unclear  that  the 
calculated  ellipse  reflects  a  related  measure  of  dispersion.  The  statistical 
properties  of  these  estimates  is  important  because  intelligence  processing 
systems  to  which  these  values  are  input  data  assume  they  are  the  estimated 
mean  and  elliptical  error  probable  from  a  bivariate  normal  distribution. 

3.6  FIX  ERROR  ELLIPSE  COMPUTATION  (CONFIDENCE  REGION) 

The  error  ellipse  computation  serves  a  dual  purpose.  First,  its 
computed  parameters  are  used  in  the  checks  to  prevent  duplications  of  previous 
fixes.  Second,  it  is  part  of  the  information  describing  a  fix. 
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The  parameters  computed  are: 


(L)  The  ellipse  orientation  angle  of  the  seraimajor  axis  in 
radians.  This  is  with  respect  to  true  north,  and  is 
defined  as  the  mean  of  all  LOB  angles  used. 

1/2 

(2)  The  ellipse  semi-major  axis  is  denoted:  a  =  s/Ctau,,)  , 

1/2 

and  semi-rainor  axis  is  denoted:  b  =  s  (tau,,)  , 

where  s  is  the  ellipse  size,  and  tau_  the  ellipse 

b 

axis  ratio  (semi-minor/semi-major ) . 


In  more  detail,  the  size  s  is  directly  proportional  to  the  root-mean 
squared  miss  angle  of  all  LOBs  used  in  fix  computation.  It  is  also  directly 
proportional  to  the  mean  (potential  weighted)  station  to  the  target  distance. 

It  is  inversely  proportional  to  a  weak  function  of  a  number  of  LOBs  used.  This 
dependence  is  and  should  be  weak  because  the  LOB  errors  in  a  ground-based 
system  are  primarily  due  to  propagation  path  perturbation  by  terrain.  These 
errors  are  not  zero-mean  and  not  uniformly  distributed. 

It  should  be  noted  that: 

(1)  The  RMS  miss  angle  is  forced  to  be  no  less  than  the 
system  instrument  accuracy  of  two  degrees.  The  actual 
value  of  the  rms  value  is  confused  in  the  available 
source  code  listings  from  2  to  8  degrees. 

(2)  The  step  size  itself  is  forced  to  be  at  least  0.5  km. 

(3)  The  step  size  is  tripled  when  the  fix  mode  is  degraded. 
This  is  based  on  the  consideration  that  with  only  two 
stations  providing  LOBs,  there  is  no  indication  of  the 
propagation  path  perturbation  effect;  hence,  a  two- 
station  fix  location  is  of  questionable  value. 
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The  ellipse  axis  ratio  semirainor/semimajor  is  defined  as  the  mean 
absolute  deviation  of  all  LOBs  used  from  the  mean  LGB  vector,  i.e.,  the 
ellipse  orientation  angle  divided  by  45  degrees. 

Note  that  when  these  error  ellipse  parameters  are  used,  they  are 
used  to  ensure  that  a  fix,  estimate  does  not  duplicate  a  p  eviously  determined 
fix.  The  criterion  employed  requires  that  the  fix  estimate  be  outside  the 
previous  fix  error  ellipse. 

3.7  COMPUTATIONAL  AND  OPERATIONAL  OBSERVATIONS 

The  TRAILBLAZER  fix  algorithm  involves  a  considerable  amount  of 
computation,  particularly  in  the  estimation  process.  The  first  fix  estimate 
requires  examination  of  every  possible  intersection  in  every  set  of  LOBs.  If 
in  the  first  estimation,  all  invalid  intersections  are  removed  from  future  con¬ 
sideration,  the  amount  of  computation  in  subsequent  estimations  (in  a  multiple 
fix  case)  is  considerably  reduced.  In  addition,  once  an  estimate  is  optimized, 
the  size  of  the  LOB  database  is  reduced  since  a  set  of  LOBs  is  used  only  once 
in  computing  a  final  fix  location.  The  fix  estimate  is  then  removed  from 
further  consideration.  Consequently,  the  amount  of  computation  in  a  multiple 
fix  case  is  reasonable. 

Since  the  entire  fix  calculation  methodology  is  highly  empirical, 
numerical  parameters  must  be  chosen  with  care.  There  are  two  k.ey  parameters: 

(1)  the  allowable  miss  angle  used  in  the  estimation  process  (currently 

3  ...igma),  and  (/.)  the  width  constant  for  the  potential  function  (currently 

2 

(1/3  sigma)  /2).  These  parameters  control  the  sensitivity  (resolution)  of 
the  fix  algorithm,  and  indirectly  affect  the  error  ellipse  size  (both  the  RMS 
miss  angle  and  mean  range  are  potential-weighted).  If  a  change  to  these  par¬ 
ameters  is  contemplated,  the  comments  in  the  code  recommend  that  their  interre- 
lationshi  -emain  reasonably  the  same.  For  example,  the  parameters  to  be 
changed  sho  .d  be  multipliers  of  3  sigma.  For  TRAILBLAZER,  the  current  param¬ 
eters  g:  aasonabie  results  according  to  the  field  test  data  (after  removing 

obviously  cad  data  attributed  to  hardware  malfunctions,  interference,  operator 
mistakes,  etc.).  These  results  indicate  that  the  error  ellipse  computed  is 
about  a  50  percent  confidence  ellipse.  A  considerable  amount  of  testing  and 
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field  experience  would  be  required  in  order  to  optimize  the  parameters,  given 
the  indeterminate  character  of  LOB  error  statistics  for  ground-based  DF.  It 
should  be  emphasized  that  this  is  a  highly  empirical  process. 

In  view  of  the  capabilities  of  the  TRAILBLAZER  system  for  LOB  set 
separation  at  collection  time  using  bin  assignments  and  subsequent  LOB 
editing,  the  fix  algorithm  should  be  used  in  the  single-fix  or  manual  mode  for 
best  results.  However,  when  no  such  separation  has  been  made  or  is  possible, 
or  when  the  operator  does  not  have  the  time  for  post-collection  LOB  editing, 
the  fix  algorithm  may  be  used  in  the  multifix  or  automatic  mode  with  a 
reasonable  expectation  of  comparable  results. 


SECTION  k 

OBSERVATIONS  AND  CONCLUSIONS 


It  appears  that  the  assembly  language  source  code  that  this  analysis 
was  based  on  might  not  be  the  latest  version.  Most  of  the  problems  indicated 
in  Section  3  with  bold-type  underlining  could  all  be  due  to  an  early  version 
of  the  software  in  transition.  The  fact  the  system  is  deployed  and  operating 
tends  to  support  the  feeling  that  the  analyzed  software  was  not  the  current 


version. 
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Rich,  Barnett.  Plane  Geometry  with  Coordinate  Geomett; 
1963. 
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1957. 

Hamming,  Richard  W.  Introduction  to  Numerical  Analysis. 
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ERROR  BUDGET 


1.0  PURPOSE 

The  purpose  of  this  appendix  is  to  identify  all  of  the  various  error 
components,  in  the  most  general  case,  when  determining  lines  of  bearings. 

These  lines  of  bearing  are  used  in  subsequent  fix  estimations  for  emitters. 

2.0  SCOPE 

The  essential  assumptions  of  this  document  are:  the  emitter  is  not 
moving  at  the  time  the  line  of  bearing  is  measured;  the  sensor  may  be  in  any 
position,  from  earthbound  to  a  moving  satellite. 

The  type  of  errors  considered  may  be  classified  into  several  cate¬ 
gories  : 


(1)  Sensor  platform  position  and  orientation  errors. 

(2)  Sensor  attitude. 

(3)  Antenna  errors. 

(4)  Instrumentation  errors. 

( 5 )  T ime . 

The  sensor  platform  postion  and  orientation  errors  may  be  referred 
to  as  "positional  errors." 

Errors  due  to  propagation  effects,  site  selection,  varying  aperture 
versus  aspect  effects,  and  operator  errors  are  not  considered  in  this  document 
Also,  errors  due  to  the  choice  of  algorithms  or  numerical  computations  are  not 
considered . 
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These  error  sources  may  be  broadly  classified  info  errors  in  the 
frame  of  reference  and  errors  in  position  measurement.  The  former  include 
errors  which  will  be  present  regardless  of  a  sensor  platform's  actual  location 
or  measurement  'hereof.  These  are  largely  the  result  of  error  or  uncertainty 
in  establishing  the  frame  of  reference  for  exchange  of  position  information. 
Position  measurement  errors  are  those  due  to  error  or  uncertainty  in  the 
methods  and  equipment  used  to  determine  platform  location  within  the  selected 
frame  of  reference. 

The  geocentric  coordinates  and  references  are: 


Latitude 

<l> 

Phi 

Longitude 

X 

Lambda 

Altitude 

h 

Orientaton  of  meridian  plane  (Direction  of  North) 

These  coordinates  are  best  described  by  the  diagram  in  Figure  1.  The  geo¬ 
graphic  latitude  is  measured  positive  from  the  equator  toward  the  North  Pole 
in  degrees.  The  geographic  longitude  is  measured  positive  from  the  prime  meri¬ 
dian  at  Greenwich  toward  the  East  in  degrees.  The  altitude  is  measured  from 
the  mean  sea  level  (the  geoid)  in  meters  and  is  positive  in  a  direction  away 
from  the  center  of  the  earth.  The  physical  sources  of  errors  in  these  para¬ 
meters  will  depend  largely  cn  the  source  of  the  data  used  to  determine  them. 
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3,1  FRAME  OF  REFERENCE  ERRORS 
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Establishment  of  a  frame  of  reference  for  exchange  of  position 
information  on  the  earth  involves  seven  major  processes,  all  of  which  may 
introduce  some  error  of  uncertainty  into  any  position  reference  within  the 
selected  framework. 
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3.1.1 


Geodetic  Errors 


Two  of  the  seven  processes  are  the  province  of  geodesy,  and  involve 
measuring  and  representing  the  shape  of  the  earth. 

One  process  involves  measurement  of  the  actual  shape  of  the  earth 
independent  of  local  variations  in  topography.  This  is  typically  most  closely 
represented  by  mean  sea  level,  i.e.,  sea  level  independent  of  variations  due  to 
lunar  tides,  and  local  gravitational  anomalies.  The  resulting  geometric 
figure,  termed  a  geoid,  becomes  the  basis  for  subsequent  representations  of  the 
earth's  surface.  This  figure  is  subject  to  error  and  uncertainty  due  to  the 
measurement  process  and  to  changes  in  the  actual  shape  of  the  earth  over  time. 

The  second  process  is  the  selection  of  a  geometrical  figure  close  to 
the  geoid  in  shape,  but  simpler  from  the  standpoint  of  mathematical  and  geo¬ 
metric  manipulation,  to  be  used  as  the  basis  of  the  mapping  process  in  a  given 
part  of  the  world.  The  figure  is  generally  based  on  a  very  nearly  spherical 
ellipsoid  which,  because  of  its  nearness  to  spherical  shape,  is  often  called  a 
spheroid.  Different  spheroids  are  in  use  for  different  parts  of  the  world  both 
for  historical  reasons  and  because  slightly  differing  ellipsoids  best  approx¬ 
imate  the  geoid  over  different  parts  of  the  earth.  Different  spheroids  are 
typically  defined  by  giving  the  radius  at  the  equator  and  the  flattening.  The 
latter  is  defined  as  the  difference  between  the  radius  at  the  equator  and  that 
at  the  pole  divided  by  the  radius  at  the  equator.  Selection  of  a  spheroid 
introduces  error  as  the  selected  figure  is  only  an  approximation  to  the  geoid, 
and  may  vary  from  the  geoid  irregularly  over  the  portion  of  the  earth  being 
mapped. 

3,1.2  Geomagnetic  Errors 

Airborne  platforms  depend  on  a  magnetic  flux  goniometer  during 
initialization  of  the  inertial  platform.  Field  soldiers  and  mobile  units  often 
have  to  depend  on  magnetic  compasses  for  determining  bearings.  Although  this 
is  one  of  the  oldest  means  of  taking  bearings,  it  can  be  very  inaccurate.  The 
earth’s  magnetic  field  tends  to  align  with  the  nearest  magnetic  pole.  However, 
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the  magnetic  poles  are  about  three  kilometers  from  the  geographic  poles. 
Furthermore,  the  two  poles,  North  and  South,  are  not  even  symmetrically  placed. 
And  to  complicate  this,  there  are  local  variations  over  all  the  earth's 
surface.  This  angle  that  the  compass  makes  with  the  grid  lines  of  a  military 
map  is  called  the  "declination"  of  the  compass.  The  magnetic  lines  of  force 
are  not  parallel  to  the  earth's  surface,  except  along  the  indefinite  circle 
called  the  magnetic  equator.  The  angle  the  magnetic  field  makes  with  a 
horizontal  plane  is  called  the  dip  angle  or  the  magnetic  inclination. 

The  declination  at  any  one  location  does  not  remain  the  same  year 
after  year  and  changes  somewhat  over  long  periods  of  time.  Besides  these 
so-called  secular  changes,  there  are  variations  within  the  year  and  also  small 
changes  of  angle  throughout  the  day.  Large  erratic  variations  occur  during 
"magnetic  storms."  These  storms  are  often  concurrent  with  the  appearance  of 
sun-spots.  Variations  from  storms  are  infrequent  enough  and  the  other  varia¬ 
tions  are  sufficiently  slow  that  it  is  practical  to  publish  maps  of  countries 
and  other  large  areas  showing  the  magnetic  declination.  On  these  maps,  points 
of  equal  magnetic  declination  are  connected  by  lines.  Each  wiggly  line  is 
labeled  with  the  amount  and  direction  of  the  magnetic  declination.  These  lines 
are  called  isogonic  lines.  The  isogonic  line  of  zero  magnetic  declination  is 
indicated  by  a  heavy  line,  and  is  called  the  agonic  line.  Maps  of  smaller  area 
indicate  the  magnetic  declination  in  their  legend  by  an  arrow  pointing  to  the 
magnetic  north  and  labeled  with  the  value  of  Che  magnetic  declination  in 
degrees . 


The  National  Space  Technology  Laboratory  at  the  Naval  Office  in  Bay 
Saint  Louis,  Mississippi  has  a  world  mathematical  model  of  the  earth's  magnetic 
field.  The  model  consists  of  an  order  12  spherical  harmonic  series  with  time 
varying  coefficients  to  take  care  of  secular  changes.  The  model  is  considered 
good  for  +  5  years.  However,  the  model  is  incapable  of  describing  anomalies 
smaller  than  about  1,100  km,  and  has  an  inaccuracy  of  about  3,500  km.  The 
model  is  updated  every  five  years  from  new  satellite  and  aircraft  survey  data. 
Local  anomalies  will  normally  deviate  a  few  degrees  of  arc  from  the  earth's 
main  field  direction,  but  can  deviate  by  tens  of  degrees  in  areas  where  the 
mineral  magnetite  is  abundant  and  in  polar  regions.  For  accurate  orientation 
using  the  earth's  magnetic  field,  there  is  no  good  substitute  for  a  local 
survey. 
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3.1.3 


Cartographic  Errors 


The  remaining  four  processes  introducing  frame  of  reference  errors 
fall  into  the  province  of  cartography,  i.e.,  the  recording,  measurement,  and 
representation  of  geographic,  topographic,  and  cultural  features  on  the  surface 
of  the  earth. 

The  first  of  these  processes  involves  selection  of  one  or  more 
coordinate  systems  to  be  used  to  specify  locations  on  the  selected  representa¬ 
tion  of  a  portion  of  the  earth.  In  virtually  all  world  reference  systems,  at 
least  one  of  the  coordinate  frames  used  will  apply  to  the  selected  spheroid, 
and  the  reference  system  used  is  in  fact  almost  always  the  familiar  geographic 
(latitude-longitude)  coordinate  system.  Errors  arise  in  this  process  due  to 
errors  in  the  measurements  associated  with  selection  of  reference  or  registra¬ 
tion  points  as  bases  of  the  coordinate  system  as  well  as  in  the  measurement  and 
computation  involved  in  extending  the  coordinate  frame  from  the  base  points 
through  the  area  to  be  mapped. 

The  second  process  involves,  in  those  cases  where  the  final  repre¬ 
sentation  will  be  planar,  a  projection  of  all  or  a  portion  of  the  selected 
spheroid  onto  a  plane  according  to  some  well  defined  set  of  mathematical  and 
geodetic  conventions.  This  step  will  often  be  followed  by  another  iteration 
of  the  first  step  to  select  a  reference  system  suitable  for  measurement  and 
computation  in  the  Euclidean  plane.  Errors  arise  in  this  process  due  to  the 
distortion  involved  in  the  projection  from  the  spheroid  to  the  plane  as  well 
as  in  any  subsequent  registration  and  extension  of  the  associated  planar 
coordinate  system. 

The  third  process  consists  of  the  recording  and  measurement  of 
surface  features  within  the  selected  coordinate  system(s).  The  errors  inherent 
to  this  process  include  those  associated  with  measurement  of  the  features  them¬ 
selves,  their  relative  locations,  and  their  locations  with  respect  to  the 
coordinate  systems  selected. 


The  fourth  process  consists  of  the  rendering  of  recorded  features 
and  associated  coordinate  systems  into  one  or  more  forms  that  can  be  inter¬ 
preted  by  people  with  a  modicum  of  training  and  experience.  Errors  arise  in 
this  process  due  to  distortions  and  simplifications  imposed  by  the  scale  and 
resolution  available  in  the  final  product,  which  in  turn  are  governed  in  part 
by  the  current  technology  and  in  part  by  the  limitations  of  the  human  percep¬ 
tual  system.  A  highway  10  meters  in  width,  for  example,  if  represented  to 
scale  on  a  1:250,000  map,  would  be  0.04  millimeters  wide  and  all  but  invisible 
to  the  naked  eye. 

3.2  POSITION  MEASUREMENT  ERRORS 

3.2.1  Inertial  Navigation 

The  four  coordinates  of  position  can  be  maintained  by  a  suitably 
designed  inertial  platform.  There  will  be  essentially  four  type  ot  errors 
with  such  systems: 

(1)  Errors  in  measurement  and  setting  of  initial  position. 

(2)  Errors  in  platform  measurement  of  inertial  change. 

(3)  Errors  in  precision  of  computation  of  position  from  inertial 
change 

(4)  Cumulative  errors  in  position,  i.e.,  drift. 

The  basic  component  of  most  modern  inertial  navigation  systems  is 
the  gyroscope.  In  addition  to  the  familiar  function  of  referencing  direction 
(gyro  compass),  gyroscopes  may  be  designed  to  measure  rotations,  to  seek  the 
local  vertical,  and  to  act  as  accelerometers. 

3.2.2  Referenced  Navigation 

Referenced  navigation  systems  are  those  that  depend  on  beacons,  or 
repeaters  ot  known  position  or  velocity.  These  may  be  classified  by  the 
geometry  of  the  data  processing: 
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The  hyperbolic  and  the  circular  navigation  systems  are  methods  of 
triangulation.  However,  the  hyperbolic  method  deals  exclusively  with  the  sides 
of  the  triangle,  while  the  circular  method  deals  with  two  sides  and  an  angle. 
The  polar  method  gives  both  a  range  and  azimuth  from  the  reference  station. 

Dacca  is  a  low  frequency  (70-130  kHz)  hyperbolic  system  that  trian¬ 
gulates  by  measuring  the  phase  difference  between  signals  from  a  master/slave 
pair  of  reference  stations.  The  master/slave  separation  is  60  to  120  n.m.  The 
useful  range  is  about  260  n.m.  over  water.  Loran  A  is  medium  frequency  (2  MHz) 
hyperbolic  system  that  triangulates  by  measuring  the  time  difference  between 
receipt  of  pulses  from  two  stations.  The  range  of  Loran  A  is  several  hundred 
miles  over  water,  but  much  reduced  over  land.  Loran  C  is  a  low  frequency 
(90-110  kHz)  version  of  Loran  A  with  considerably  more  range.  OMEGA  is  a  very 
low  frequency  hyperbolic  system  that  triangulates  by  comparing  the  phase  of 
signals  from  two  beacons  separated  by  a  baseline  of  5,000  to  6,000  miles.  The 
coverage  is  world  wide  and  may  be  used  by  submersibles . 


II 
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Satellite-aided  navigation  has  the  most  diverse  possibility  for  use 
as  a  referenced  system  of  navigation.  The  orbital  ’lements  and  thus  both  the 
position  and  velocity  of  the  satellite  are  accurately  known.  By  combining  such 
measurables  as  elevation  angle,  azimuth  angle,  ranges,  difference  in  range, 
range  sum,  or  doppler  shift,  fixes  may  be  obtained  that  fit  any  of  the  listed 
categories  in  the  first  paragraph  of  this  section.  Methods  that  depend  on 
measurement  of  the  elevation  angle  of  one  or  more  satellites  determine  small 
circles  on  the  earth's  surface  for  fixes.  Methods  that  determine  distances 
lead  to  hyperbolic  conic  lattices  for  fixing. 


TACAN  is  a  UHF  radio  navigation  system  which  provides  both  distance 
and  bearing  information  of  the  aircraft  relative  to  the  selected  ground  beacon. 
The  antenna  system  is  the  key  to  measuring  the  aximuth.  The  antenna  system  has 
a  single,  central  element  for  transmission  and  reception.  Tie  parasitic 
elements  are  mounted  on  two  concentric  cylinders  which  rotate  at  fifteen 
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revolutions  per  second.  The  inner  cylinder  consists  of  a  single  parasitic 
element  which  causes  a  single  cardioid  polar  pattern  rotating  at  15  rps.  The 
outer  cylinder  has  nine  parasitic  elements  that  superimpose  nine  lobs  on  the 

cardioid  pattern.  This  gives  an  amplitude  modulated  signal  with  two  frequency 

modulations  of  15  Hz  and  135  Hz.  The  transponder  further  emits  bearing  refer¬ 
ence  pulses  as  the  peak  of  each  lobe  points  East.  When  the  lobe  which  coin¬ 
cides  with  the  peak  of  the  cardioid  points  East,  a  special  "North”  reference 
pulse  code  is  transmitted.  The  airborne  equipment  measures  the  phase  relation¬ 
ship  of  the  maximum  signal  amplitude  relative  to  the  North  reference  pulses  in 
order  to  determine  the  bearing  of  the  aircraft  relative  to  the  beacon.  The 
accuracy  of  the  azimuth  is  in  the  order  of  magnitude  of  two  degrees.  The 

distance  measuring  part  of  TACAN  equipment  is  like  radar  except  that  the  return 

signal  comes  from  a  beacon  used  to  produce  strong  artificial  echoes.  The 
beacon  will  respond  to  numerous  simultaneous  interrogations.  To  make  this  pos¬ 
sible,  the  pulse  repetition  rate  of  the  airborne  transmitter  is  cause  to  jitter 
in  a  random  manner.  The  receiver  is  allowed  to  recognize  only  those  pulses 
received  that  follow  the  same  jitter  pattern  and  ignore  all  other.  The  slant 
range  is  determined  to  roughly  0.25  nautical  miles  under  most  conditions. 

3.2.3  Doppler  Navigation 

Airborne  Doppler  is  a  SHF  (micro-wave)  system  of  navigation  using 
the  terrain  or  water  below  as  a  reference.  Depending  on  the  particular  doppler 
system  used,  some  or  all  of  the  following  data  can  be  made  available  to  the 
crew: 

(1)  Component  velocities  and  distances  run,  along,  across,  and 
perpendicular  to  the  aircraft  axes. 

(2)  Ground  speed. 

(3)  Drift  angle. 

(4)  Angle  of  attack. 

(5)  Height  above  terrain. 


If  True  Air  Speed,  Pitch,  and  Heading  Angles  are  available  from 
such  sources  as  an  inertial  system,  then  the  following  secondary  data  may  be 
obtained. 

(1)  Wind  speed  and  direction. 

(2)  Climb  angle. 

(3)  Track  angle. 

The  Doppler  systems  may  have  various  configurations  of  antenna  beams 
directed  at  different  angles  toward  the  earth.  A  two  beam  system  may  be  used 
to  measure  ground  speed  and  drift.  A  three  beam  system  is  basically  sufficient 
to  extract  ail  three  velocity  components,  but  a  four  beam  symmetrical  arrange¬ 
ment  is  often  used. 

4.0  ATTITUDE  ERRORS 

4.1  ATTITUDE  COORDINATES 

The  three  attitude  coordiuaces  are; 

Roll  Angle  a  alpha 

Pitch  Angle  B  beta 

Yaw  Angle  Y  gamma 

Figure  2  serves  to  define  each  of  these  angles.  These  are  the 
standard  Euler  angles  as  defined  by  a  "right  hand"  rule.  However,  it  should 
be  noted  that  the  sign  of  these  angles  vary  considerably  throughout  published 
literature.  See  Korn  and  Korn,  reference  2,  section  14.L0-6,  for  a  discussion 
of  this  coordinate  system  and  the  diverse  choice  of  signs.  In  some  aireborne 
systems  these  positional  coordinates  are  limited  by  preset  stops  which  may 
introduce  non-linear  errors. 
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With  a  cursory  examination  of  these  six  coordinates,  it  is  apparent 
that  errors  in  three  of  them  will  produce  the  larger  errors.  An  error  in  yaw 
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angle  alone  will  produce  a  divergence  of  the  azimuth  angle  of  bearing.  This 
azimuth  error  will  always  be  quite  close  in  magnitude,  but  opposite  in  sign,  to 
the  Yaw  error.  Errors  in  longitude  and  latitude  will  produce  an  error  in  the 
position  of  the  line  of  bearing  «s  a  function  of  the  azimuth  angle,  but  this 
does  not  effect  the  azimuth  angle.  If  the  azimuth  angle  is  in  the  vicinity  of 
zero  or  180  degrees,  an  error  of  longitude  will  be  reflected  directly,  and  of 
nearly  the  same  value,  in  the  longitude  of  the  fix  estimations.  At  azimuth 
angles  of  90  and  270  degrees,  the  line  of  bearing  and  consequently  the  fixes 
are  unaffected  by  errors  in  longitude.  The  effects  of  errors  in  latitude  are 
analogous  in  their  effect  but  displaced  by  90  degrees. 

It  is  not  so  obvious  that  an  error  in  the  three  remaining  coor¬ 
dinates  (altitude,  roll,  and  pitch)  should  have  any  effect  on  the  line  of 
bearing.  Indeed  an  error  in  altitude  alone  should  only  change  the  slant  range 
and  have  no  effect  on  the  line  of  bearing.  However,  when  coupled  wih  errors  in 
roll  and  pitch,  there  is  a  definite  mathematical  relation  or  coupling.  The 
significance  of  error  in  altitude  remains  to  be  evaluated.  Errors  in  roll 
and  pitch  (which  have  less  effect  on  the  error  of  the  fix  estimate  than  yaw, 
longitude,  and  latitude)  directly  cause  errors  in  azimuth  angle  on  the  line  of 
bearing . 


5.0  ANTENNA  ERRORS 

Orientation  with  respect  to  the  platform. 

Difference  between  the  mechanical  axis  and  RF  axis. 

Beam  width. 

These  first  two  antenna  errors  are  directly  related  to  the  platform 
attitude  coordinate  errors.  In  fact  the  orientation  of  the  antenna  with 
respect  to  the  platform  and  the  difference  between  the  mechanical  axis  and  the 
RF  axis  are  best  described  by  Euler  angles.  If  the  axis  defining  these  coor¬ 
dinates  are  chosen  originally  in  coincidence,  first  order  approximations  will 
sptvp  tn  considerably  simplify  the  maze  of  t ri gonome trie  functions  relating 
these  angles.  These  three  Euler  angles  can  be  identified  as  pitch,  roll,  and 
yaw.  For  small  errors  in  these  angles,  the  errors  may  be  simply  added  to  the 
corresponding  platform  angles.  It  should  be  noted  at  this  point  that  the  RF 
boresight  error  is  a  function  of  the  radio  frequency. 
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Beam  width  is  always  a  function  of  the  antenna  geometry  and  fre¬ 
quency.  A  phased  antenna  system's  beam  width  will  vary  considerably  with 
change  in  aspect  angle. 

6.0  INSTRUMENTATION  ERRORS 

Bias  (Systematic  errors). 

Noise  (Random  errors). 

Bias  errors,  for  example,  are  systematic  errors  such  as  boresite 
errors,  parallax  errors  in  instrument  readings,  and  bezel  errors.  Bias  errors 
are  usually  minimized  by  calibration  procedures. 

Noise  errors  are  due  to  random  phenomenon  such  as  receiver  noise. 
This  noise  normally  produces  random  errors  in  bearings  by  increasing  the  region 
of  uncertainty  when  determining  the  minima  of  a  signal  or  the  change  in  sign 
from  the  phase  of  a  signal.  There  are  many  techniques  of  minimizing  the 
effects  of  noise,  depending  on  the  source  and  nature  of  the  noise  (see  Refer¬ 
ence  5).  In  high  frequency  receivers,  the  receiver's  "front  end"  is  a  high 
source  of  thermal  noise,  so  the  high  gain  required  is  usually  obtained  after 
heterodyning  to  a  lower  frequency  or  after  further  detection  at  the  "rear  end." 
Commonly,  the  band  pass  of  filtering  is  reduced  to  the  minimum  that  will  not 
deteriorate  the  informatio'  content.  The  effect  of  impulse  noise,  such  as 
noise  emanating  from  electrical  ignitions,  can  be  minimized  by  amplitude 
clipping  just  above  the  signal  level. 

7.0  ERROR  TABULATION 

The  sensor  positional  error  is  equally  important  in  fixing,  mobile, 
or  airborne  sensors.  The  attitude  errors  are  most  important  in  airborne 
sensors.  The  sensor  geometric  error  refers  to  such  errors  as  the  difference 
between  the  geometric  and  RF  axis  of  a  direction-finding  antenna,  or  even  an 
optical  tracker.  Range  is  included  with  the  geometric  sensor  errors  for  con¬ 
venience  only. 
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The  specifications  and  tolerances  will  always  include  the  units. 
The  exact  meaning  of  the  specification  and  tolerance  columns  will  depend  on 
the  instrument  involved. 
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COMPONENT^ 
A/N  NUMBER 
MODEL _ 


CLASSIFICATION 
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1.  Sensor  Positional  Errors 

a.  Longitude _ 

b.  Latitude _ 

c.  Altitude _ 

d.  Position  (linear  distance ) 

2.  Sensor  Attitude  Errors 

a.  Reference  meridian  (North) _ 

b.  Roll _ 

c.  Pitch _ 

d .  Yaw _ 

e.  Rates  (TBD) _ 

3.  Sensor  Geometric  Errors 

a.  Azimuth _ 

b.  Elevation _ _ _ 

c .  Range _ _ 

4.  Instrument  errors 

a.  Bias  (systematic  or  secular  errors) 

b.  Noise  random  errors _ 
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Figure  2.  Attitude  of  the  Platform 
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APPENDIX  C 

TRAILBLAZER  ERROR  BUDGET 


TABULATION  OF  ERROR  BUDGET 


SYSTEM _ 

COMPONENT 
A/N  NUMBER 
MODEL 


CLASSIFICATION 


1.  Sensor  Positional  Errors 

a.  Longitude _ 

b.  Lat itude _ 

c.  Altitude _ __ 

d.  Position  (linear  distance) _ 

2.  Sensor  Attitude  Errors 

a.  Reference  meridian  (North) _ 

b.  Roll _ 

c.  Pitch _ 

d .  Yaw _ 

e.  Rates  (TED) _ 

3.  Sensor  Geometric  Errors 

a.  Azimuth _ 

b.  Elevation _ _ ___ 

c .  Range _ 

A.  Instrument  errors 

a.  Bias  (systematic  or  secular  errors) 

b.  Noise  random  errors 


MEASURED  LOBs 


PROGRAM  TRAILBLAZER(  INPUT*  OUTPUT*  TRAILIN.  TRAILOUT); 

********  ****************  ******** ************  ********************* 

<* 

Thi»  program/algor ithm  has  been  tailorad  specifically  for  operation 
with  tha  LOB  data  produced  by  a  ground— based  DF  NET*  consisting  of  a 
small  fixed  number  of  DF  stations,  whose  locations  remain  invariant 
during  the  data  collection  process..  The  distinctive  feature  of  this 
algorithm  is  its  capacity  to  resolve  reliably,  multiple  targets 
while  at  the  same  time  avoiding  GHOSTS,  i.  e.  false  targets  arising 
from  coincidental  intersections  of  unrelated  LOBS. 

Original  program  written  in  " ROLM  ASSEMBLEY  LANGUAGE". 

Translated  into  PASCAL  by  Nick  Covella.  JUNE  1984. 


CONST 


tor  =  3* 


Shap  efac 
Pie 

Radian 

Max 

Min 

Sigma 

TYPE 


■  3.  14159; 
»  180; 

«  101; 

*  0.  6; 

=  8; 


Stati ons 


Data 


=  RECORD 


Xcoord 

REAL;  ( * 

X-Location  of  any  station.  *) 

Ycoord 

REAL;  (* 

Y-Location  of  any  station.*) 

Theta 

REAL;  ( * 

Angle  inputted  from  the 
user.  *) 

Alpha 

REAL;  <  * 

Calculated  angle  from  input 
to  TRUE  NORTH  *> 

MissAng 1 e 

REAL;  ( * 

Difference  between  the  Alpha  angle 
and  the  recalculated  angle  after 
an  intersection  has  been  found.  •«- ; 

Semimajor 

REAL;  <* 

in  kilometers  *) 

Semiminor 

REAL;  <* 

in  kilometers  *) 

Di stanc  e 

REAL;  <* 

approx,  distance  from  the  staion; 
to  the  object  being  "fixed".  *) 

Ref  1 ob 

BOOLEAN; 

(*  indicates  that  an  LOB  has 

been  attempted  by  the  station. 

Flag 

BOOLEAN; 

<*  indicates  that  an  intersection 
for  this  station  has  been 
calculated.  *) 

Org  i  n 

END; 

INTEGER; 

(*  indicator  for  station 
manipulation  *) 

-  RECORD 

Xintercept 

:  REAL; 

Yintercept 

:  REAL; 

Support  d 

—4  REAL; 

Sup  o  ortF 1 ag 

:  ARRAY  C 1 

,  53  OF  BOOLEAN; 

SupC  ountF  lag  ;  INTEGER;  _ 

Flag 


BOOLEAN; 


END; 

Data  structure  that  will 
keep  track  of  the  number 
of  LOB ' s  for  the  station 
being  accessed.  *■) 


*  Table  for  the  intersections 
of  the  lobs  in  each  bin.  * ) 


VAR 


Missed  Ang 1 e 

REAL, 

AS  t  at l on 

INTEGER 

13  i  nNumb  er 

INTEGER 

Index 

I NTEGER 

Intlosuhscript 

I NTEGER 

I n  tPas  s 

INTEGER 

Nma  x 

INTEGER 

Nether 

I NTEGER 

Ns  ame 

INTEGER 

Col 

INTEGER 

S  t aNumb 

INTEGER 

Row 

INTEGER 

TrailQut 

TEXT, 

Tr  a  i  1 1  n 

TEXT, 

I  n  t  lo 

I ntlotype 

Lob 

Lobtu  pe. 

ID 

Stat  i  ons; 

Table 

TVT, 

Intlotype  =  ARRAYE1.  53  OF  INTEGER;  (* 


Lobty pe 

=  ARRAYC1. 

.  53 

OF  Stations; 

TVT 

=  ARRAYC  1. 

.  5. 

1.  .  53  OF  Data;  ( 

r  *** *****♦■»******■*******•»«-***•*-»**■«-»* *♦*■»*****■**«•  *****^*****  *****#■*■»• ■*■**■*■*■*«■■*  v 
PROCEDURE  CLEAR ( VAR  Lobrec  :  Lob  type; 

1  I  AD  T  1  «  i..  ^  n  a- 

v  r-n  \  i,n  v  i  uai  i  ai)  ill  v  iu 

VAR  Tabletype  TVT  ) ; 

<  **#  *■*-**  ****■■#*#■*  *^**-IH»*******HI-#*^*******a«-#**a»**#**-IHHt'*  ■»*##  **•■#*  *■*■**  ***•*■ 

<*  This  subroutine  clears  all  of  the  data  structures  and  prepares  them 
for  either  the  first  pass  or  any  other  subsequent  passes. 

*  i 

VAR 

J.  INTEGER; 

BinNumber  INTEGER; 

AS  tat i on  INTEGER; 

I  n  1 1  os  ub  sc  r  i  p  t  INTEGER; 

3  EG  I  N 


fi 


WRITELNC ' 

I  =0, 

FOP.  B  l  nNum b  er 


BEGIN 


Entering  CLEAR'); 

TQ  5  DO  <*  initialize  the  variables  of  the  bins  of 
each  station.  *  > 

D-5 


Lo  brec 
Lobrec 
Lobrec 
Lo  brec 
Lo  brec 
Lo  brec 
Lobrec 
Lobrec 
Lobrec 
Lobrec 
Lobrec 


CB  1 n  Numb 
CBinNumb 
CBinNumb 
CBinNumb 
CBinNumb 
CBinNumb 
CBinNumb 
CBinNumb 
CBinNumb 
CBinNumb 
CBinNumb 


Xcoord 
Ycoord 
Theta 
Alpha 
MissAng 1 e 
Semi  major 
Semi  minor 
Dist  anc e 
Ref  1 ob 
Flag 
Org  i  n 


=  00.  000; 
=  00.  000; 
=>  00.  000; 
=«  00.  000; 
*  00.  000; 
*  OO.  000; 
=  00.  000; 
=  00.  000; 
»  FALSE; 

»  FALSE; 

=  0; 


FOR  BinNumber  : =  1  TO  5  DO 


BEGIN 


I  :  =  I  +  1, 

FOR  AS  tat  ion  ; 


=  1  TO  5  DO  <*  initialize  the  truth-table  of  the  bins  + 


BEGIN 


Table typeCBinN umber. AStationl.  Xintercept 
Ta b let y peC B inNumb e r , AS tat i on 3 .  Vintercept 
Ta b  1  e t  y  pe  C  EinNumber. AS tat i on! .  Support 
TabletypeC  B inNumb  e  r . AS  tationl.  SupCountFlag 
TabletypeCBinNumber.AStationl.  SupportFlagC II 
TabletypeC  B inNumbe  r , AStationl,  Flag 


=  00. 000 
=  OC .  000 
=  00 . 000 
=  Oi 

=  FALSE; 
=  FALSE; 


FOR  Intlosubscnpt  ;  =  1  TO  5  DO  (*  intersection  LOB  array  *■) 

I  nt  1  oar  r  ay  C  I  n  tl  o  sub  s  cr  ip  1 1  :  =  0;  (*  initialized  to  -1  in  actual  progi 


WR I TELN ( ' 


Leaving  CLEAR'); 


i*  PROCEDURE  "CLEAR'1  *) 

■*****#*»  #»*#»»***•»#***######  ♦♦*#*##*#####***#***-»*-t**#^4*  +■■*#■***-*  * 


:EDURE  I  NELL  I  P  SE  <  VAR  Xefix 

VAR  Yef l x 
VAR  Answer 
Table 


REAL; 

REAL; 

BOOLEAN; 

TVT  > ; 


tt-**-*-******^*  ■*■«-***■*-*****  ■*****#**■*****#■*-**■*■*-*■**■«••*■*•*■«•-»•*  •*■**■*•  +  *-■ 

his  procedure  will  insure  that  a  fix  estimate  doesn't  duplictate 
:he  last  or  any  previous  fix  INELLIPSE  returns  true  in  th  e  b  oolean 
ariahie  Answer  if  the  n  x  estimate  is  in  the  ellipse  determined  oy 
;he  previous  fixes. 


Va  )  ue  1  .  Value iSqr,  Va  1  ue2. 
Tempaddi  Tempcos,  Ternpsin 
Sum 

■d _ 


Va  1  u  e2S  qr 


REAL; 

REAL. 

REAL, 


REAL; 

REAL; 


V.  S 9. tu mi  n ora  x  is 

v 

ijenumajoraxis 


BEGIN 


WR  I  TELN  i  ' 


Entering  INELLIPSE'); 


<  *  get 
( *  get 
Tempco 
Tempsi 
Temp ad 
Semima 
Valuel 
Va  1  ue  1 
T  e  m  o  c  e 


jora 


Temps:  i 
T  e  m  n  a  d  - 


Sem imi 
Va 1 ue2 
Va  I  ue2 
Sum  :  * 
IF  Sum 
Answ 


Sqr 
Va  1 


nd  y  coordinates  of  the  center  of  the  ellipse  *) 
cosine  of  the  ellipse  *) 

COS ( ( Ye f i x  *  Pie)  /  Radian); 

SIN  < (Xefix  *  Pie)  /  Radian); 

Tempcos  +  Tempsin; 

xis  :  =*  Tempcos  +  Tempsin*  <  *■  JUST  FOR  ARGUEMENT  *) 

Temp  add  /  Semima jora x is; 

:  =  Value!  *  Value!. 

COS ((Xefix  *  Pie)  /  Radian), 

SIN(CYefix  *  Pie)  /  Radian); 

Tempcos  —  Tempsin; 

xis  ; =  Tempcos  ■+•  Tempsin;  I#  JUST  FOR  ARQUEMENT  * ) 

Tempadd  /  Semimi nora x is; 

:=  Value2  *  Value2; 
uelSqr  +  ValueSSqr; 

THEN 

=  TRUE  (*  indicating  that  the  estimate  mas  probably  d 
from  any  other  estimate  *) 


<*  JUST  FOR  ARGUEMENT  *) 


i  f  f  e  r  e  n  fc 


ELSE 


•  d  •  X  »  —  L  L  ^  L  L  U  • 

i  u  l  u  ci  v*  a  il  y  w  1 1  a  v  V  li  x 


WR  I  TELN  ( 


Leaving  INELLIPSE'); 


.-■END;  (*  PROCEDURE  “I  NELL  IPSE"  *) 


N 


PROCEDURE  MS CAN ( VAR  Lobrec  :  Lob  type); 

\  -rt- ■!<*■«* -*!■  ’JHHHt ♦  ■*##•&  -ft- -fr 

<  * 

This  procedure  searches  all  LOBS  in  the  local  database,  if  the 
LOGS  have  been  marked  indicating  that  they  have  been  used  for  a 
particular  fix  then  they  mill  be  unmarked  by  this  procedure. 

*> 


r-  -&■  -x  •>  r-  >  - . 


BinNumber  INTEGER; 


;EGI  N 


WR I TELN  < ' 


Entering  MSC^J'); 


FOP  BinNijmber  :=  1  TO  5  DO  (*  check  each  BINSET  to  see  which  have  been 

marked.  Unmark  those  that  have  been  markea 


IF  Lobr  ec  CB  mNumber  3  FI  ag  =  TRUE  THEN 
LobrecCBmNumber],  Flag  :=  FAt_SE» 


"  WR I  TEL  N  ( 

.a 


Leaving  NS CAN ' ) , 


END •  (*  PROCEDURE  "MSCAN"  *) 


PROCEDURE  SCNDB1'  Cornu  uit  INTEGER, 


i 

/v 

<%■ 

r. 

■f. 

i 


f? 


$ 


I 

s 


Reflob  :  INTEGER; 

VAR  Lobrec  Lob  type), 

(  *********************************** ***##«***♦*##*♦***■»#**#***#♦*+**  it  *#**■»+■*•» 
<  * 

This  procedure  accumlates  counts  of  marked  LOBS  for  eoch  station  in 
the  counter  NOTHER  and  NSAME.  "SCNDB"  is  called  by  "FPQT"  to  per¬ 
form  miS5angle  variance  on  a  station  basis. 

* ) 


£ 

>■  ;■ 

l  Vi 


VAR 


BEGIN 

WR  I TELN  ( 


Entering  SC NDB ' ) ; 


a  :  =  C  omp  c  n  t, 

Y  .  »  Ref  lob, 

SumVariance  :=  Q; 

SumPotential  O; 

No  th  er  : =  O; 

Nsame  =0; 

IF  X  =  1  THEN 

Scancntr  =  Nother; 

For  BinNumber  : =  1  TO  5  DO 
IF  Lobret CBinNumberl.  Reflob  =  FALSE  THEN 

BEGIN 

Lobr  e  c  IB  i  nNumber  ]  .  Fla  g  :=  FALSE; 
Scancntr  .=  Nsame. 

Nother  : =  1 ; 

END . 

Lo b rec C B inNumb er 3  .  Reflob  .  =  TRUE; 

WHILE  Lobr ec CB inNumber 3  Flag  O  TRUE  DO 
BEGIN 

IF  X  =  1  THEN 

Scancntr  :=  Scancntr  +  1, 

IF  ((Scancntr  =0)  OR  (X  =  0))  THEN 
Hissed Angle  :=  Scanc  ntr  -  X  (*  LOB 
IF  HissedAngle  '=  O  THEN 

Hissed Angle  : =  -(HissedAngle); 

IF  HissedAngle  >  Pie  T HEN  n  Q 

D~o 


-  reflob  * )  ; 


.V 

<3 


rs 

L. 


SumV'ar  iance 

REAL, 

Su  mPo  t  ent i a  1 

REAL; 

Tempi.  Temp2-  Temp3 

REAL, 

»i 

1  ^ 

Ei  p 

REAL; 

V- 

r  o  t  en  t  i  a  1 

REAL, 

B l nNum  b  er 

INTEGER; 

Sigma.  Shapefacfcor 

INTEGER; 

(* 

constants 

given  in  the  progeam  *) 

£ 

X,  Y 

INTEGER; 

<* 

Sub  s  t itute 

variables  for  the  parameters 

passed  to 

this  subroutine  *> 

Sc  anc  ntr 

INTEGER; 

<* 

counter  of 

marked  LOBS  *) 

k 

N" 


R 


ft 

bv 


*  * 
I  a 

l  - 


BEGIN 


Mi ssedAngle  (2  *  Pi*)  -  MissedAngle; 

Sigma  : ■  8; 

Shapef actor  :  ■  3; 

T*mp 1  Sigma  »  Shapefactor; 

Ttmp2  :»  MissedAngle  /  Tempi; 

Exp  :»  (  (SQR  (T*mp2) )  /  2); 

Temp3  : ■  (((Exp  *  Exp)/2)  +  ((Exp  *  Exp  *  Exp>/6>  +  Exp  +  1); 
Potential  :■  1/Temp3» 

SomPotantial  :■  Sumpotantial  +  Potantial; 

SumVariance  :■  SumVariance  +  ( (SQR (MissedAngle )  *  Potential)); 

END; 

END. 

IF  Nother  =  0  THEN 
Scancntr  : =  Nother 
ELSE 

Scancntr  :=  Nsame; 

WR  ITELNC  '  Laaving  SCANDB  ' ); 

D;  (*  PROCEDURE  "SCANDB"  *) 

»•****•* *#**•*•**•***•«•* *■*■»•***** *****HHH»#*******#*****#**********tf**-*i****-S**  £■»>•■*■>" 


OCEDURE  XPREP <  Stldent  :  Lobtype; 

I  :  INTEGER; 

Int  x  :  REAL; 

Inty  :  REAL; 

VAR  Computed  :  BOOLEAN); 


****««*-»«^  *«***«*«  ***«*«-»*  *  *^###a.*******e  •**■**•**■**■*  ***•»•***■*  ************-*■**■-**■>' 

This  procedure  computes  the  LOB  angle  of  the  Intersection  and 
returns  the  coordinates  of  the  proper  location  in  the  variables 
Intx  and  Inty. 


4* 

DistanceX  .  REAL; 

OistanceY  REAL; 

TotalDist  ;  REAL; 

:qin 

WRITELNf'  Entering  XPREP'); 


DistanceX  ■*  S  t  Id  en  tC  1 1 .  Xc  oor  d  -  Intx; 

DistanceY  ■  St  IdentC  I  3.  Yc  oord  -  Inty; 

To  talD x  st  -  SQRT( (SQR (Distanc *X> )  ♦  ( SQR ( Di s t anc e Y ) ) )  ; 

IF  (TotalDist  C  Min)  THEN 
DEC  I  N 

WR  ITELN'.  Trai  lOut)  .  D"9 

WP  I  TELN '  T  -  a  i  lOu  t.  'The  distance  from  Station#'.  St  Id  entCIl.  Or  gin:  1.  i?  tc 
WR  I  TELNf  T  •*  a  \  l  Out.  close  to  obtain  a  proper  intersection.  '); 


WR I TELN (TraxlOut); 

In  t  x  :  =  00.  000;  ! 

Inty  :  =  00.  COO; 

Computed  . =  FALSE.  t 

END; 

IF  <To talD i st  >  Max)  THEN  ’ 

BEGIN 

WR IT£LN(Tr ai lOut ) ;  j 

WR  I  TEEN  (TrailOut*  'The  distance  from  Station  #  ' ,  St  I  d  ent  C  I  3  .  Orgin  i-  ii  ■- 
WR  I  TELN  (TrailGut.  'far  to  obtain  a  proper  intersection.  ');  j 

WR  ITELN ( Tr a l lOu t ) ;  j 

In  t  x  :  =  00.  000; 

Inty  : =  00 . 000 ; 

Computed  ; =  FALSE;  1 

END; 

WR I TELN ( '  Leaving  XPREF ' ) ; 


END;  <*  PROCEDURE  "XPREP"  *) 


(  *•**■**■#***■*■#•*  V*^**#**********^******^*  a-*********"*-*-*#**-*-*-**-*  *-*-*-*-<•*'*  <<- 


PRuCEuURh  XCunPv 


I 


I N  i  tGER; 

INTEGER; 

REAL; 

REAL; 

Lobty  pe; 

REAL; 

REAL; 

BOOLEAN) ; 

(  «■  This  procedure  is  called  twice  by  "XSEC"/  once  for  each  of  the  two 

in  the  computation  of  their  LOB  intersection.  If 
then  no  intersection  is  computed. 


VAR 

VAR 


Xdist 
Ydist 
St 1 d ent 
Int  s 
Inty 


VAR  Computed 


stations  involved 
the  intersection  is  negative 


NOTE.  The  parameters 
stations  whose 


* ) 


VAR 


of  I  and  J  contain  the  staion  number  of  the  two 
intersection  is  to  be  computed. 


BEGIN 


Dy 

REAL;  <*  Temporary  variable  for 

in  e I .  Sine  I 

:  PEAL; 

ineJi  SineJ 

:  REAL; 

pD  i  s  t 

:  REAL; 

TELN  (  ' 

Entering  XCOMR ' ) ; 

--  Xdist. 

:  =  Ydist; 

(St  IdentCII. 

Alpha  =  StI  dent  I.  J  j.  Alpha)  THEN 

i£kJ 


BEGIN 


D-10 


Intx  :  =  00 .  000 ; 

Inty  :  =  00 .  000 ; 

Computed  :  =  FALSE; 

WR.  ITELNCTra  i  lQut ) ; 

WRITE(Trai  1  Out,  'Station  #  '  ,  St  I  dent  C 1 3  Org  i  n ;  1 ,  '  and  Station  #').• 
WR  ITE (  Tra  i lOut,  StIdentCJ3.  Org in: 1 ,  '  have  the  same  initial  LOB'); 
WR  ITELN( Trail Out) ; 

WR  ITELNCTra  ilOut,  '  traj  ectory.  ...  No  intersection  possible.'); 
WRITELN(TrailOut); 


ELSE 


BEGIN 


Co  sine  I 
CosmeJ 
S i ne I  : = 
Si  ne  J  ; = 
Xdist  :■ 
Xd  i st  :  = 
Yd ist  : = 
Yd i st  ; = 
Computed 


=  COS  (  ( <  St  I  d  en  tC  I  1 .  Th  eta  )  *■  Pie) 
=  COS  (  ( (Stldentc  Jl.  Theta)  *  Pie) 
SIN( ( (StldentCU.  Theta)  *  Pie)  / 
SIN< ( (StI dent CJ3.  Theta )  *  Pie)  / 
<  (  Dx  *  CosineJ  )  —  (  Dy  *  SineJ 


/  R  a  d  i  a  n  ) ; 
/  Radian), 
Radian); 
Rad i an ) ; 

)  ); 


Xd  ist/SIN(  ( (StIdentCJ3.  Theta  -  StldentCII. 
((  Dx  *  Cosinel  )  -  (  Dy  *  Sinel  )); 

Yd  ist/SIN< ( (StldentCJl. Theta  -  StIdentCI3 
:  =  True, 


Th  e  ta ) 


Theta) 


Pie) 


Pie) 


Rad  ia-- 


Rad  is- 


IF  (Computed  =  TRUE)  THEN 


BEGIN 

Intx  :=  St  I  den t C 1 3 .  Xc o nrd  +  (Xdist  *  SIN( < ( ( St  I  den t [ I  3 . Th e t a > 

*  Pie)  /  Radian))); 

Inty  ;=  St  I  den t C I  3 .  Yc o ord  +  (Xdist  *  COS( ( ( (St  I  den t Cl  3 . The ta ) 

*  Pie)  /  Radian)))i 
Computed  ;=  TRUE; 


IF  ((ABS(  (StldentC  13.  Alpha) 
(Computed  ~  TRUE) )  THEN 


(St  I  dent  CJ3.  Alpha)  >  O 


BEGIN 


WR I TELN (TrailOut) , 

WR  ITE( Tra i lOut ,  'The  difference  between  the  angles  of  Station  #'); 

WR  ITE ( Tra i lOu  t , St  I  dent  C 1 3 .  Org i n ;  1 ,  '  and  Station  # ' .  St  I d  ent  C J3 .  Or g i n :  1 )  ; 

WR I TELN (TrailOut), 

WR ITELN ( Tr a i lOu t,  'is  too  small  to  obtain  a  proper  intersection.  '), 

WR  ITELN(Tr  a  ilOut,  'Station  #  ' ,  S 1 1  d  e  n  1 1 1  3  ,  Or  a  in  ;  1,  '  and  Station  #',StIdent! 
WR ITELN( Tr a i lOut,  'are  not  participating  in  the  test  data'), 

WR  I  TELN (TrailOut)  i 
Intx  : =  00  000; 

Inty  ■  =  00  000 ; 

Computed  :  =  FALSE, 


BEGIN 

XPREP  (Stldent,  StldentCIl  Qrgini  In  t  x ,  Inty.C  omp  u  ted  )  ; 
XPREP(StIdent>  StldentCJI,  Orgin,  In  t  x ,  Inty>  Computed  )  i 

END; 

WR  I TELN (  '  Leaving  XCOMP'); 

END;  (*  PROCEDURE  XCOMP  *) 


(  tt*###********-**  tt-B"**  ##*#*«*#  *****#**#**■£■«••#•*•#■*■«■##  ******•*■«■  *-H'«  ***<r 


PROCEDURE  I LOSS  (  VAR 

AllStat 

Lob  type;  <  •*  LQB 's 

back  to  remaining  stati 

which 

have  inputted  data  into 

the  da  taba  se.  -#■) 

Stldent 

Lo  b  ty  p  e; 

I  lit  X 

REAL;  ( *  possibly 

a  VAR  *) 

Inty 

REAL;  <  *  possibly 

a  VAR  *) 

StaNumb 

INTEGER; 

VAR 

Table 

TVT); 

(  *•**  •«■**■» •«•■&•#•*■  #***  •***■*  ******■»*•«■•#■*#***■*■#•■*#■■*#■#•**#*■**■«■-#■#-*•«■#■•**  *•**■* #•**-**■*  * -ft-  is-**.# 

(*  This  procedure  completes  the  intersection  file  called  "INTLO". 

Reported  LOBS  are  flagged  and  stored  in  "INTLO"  as  indexed  bg  the 
repec  tive  station  numbers. 

■*  i 

VAR 


An  g leTheta 

REAL; 

(  * 

Angle  formed  by  the  fixed  point  ba 
station.  * ) 

c  k  to  ea h 

Di  stX 

REAL, 

( * 

distance  between  the  fixed  x-coord 

and  the 

x-coord  of  the  stations'  position. 

*) 

Di  stY 

REAL, 

(  * 

distance  between  the  fixed  y-coord 

and  the 

g-coord  of  the  stations'  position. 

*) 

Te  mp  A 

REAL; 

Te  mp  5 

REAL; 

I.  J,  K,  n 

INTEGER; 

A,  0 

INTEGER; 

Su  p  p  or  tCoun t 

INTEGER; 

Ne  wDa  t  a 

Lob  type; 

BEGIN 


WR I TELN ( '  Entering  ILOBS'); 

N  :  =  1  , 

Di  stX  :  =  00.  000; 

Di  stY  :  *  00.  000, 

K  : =  S  taNumbi 

FOR  J  :=  1  TO  StaNumb  DO 

IF  <  (Stldent  CJ3.  Ref  1  ob  =  TRUE)  AND  (StldentCUJ.  Flag  =  TRUE))  THEN 
BEGIN 

New Data  CN3  .  Orgin  :  =  St  I  dei.t  C  Jl.  Qt  g  in; 

N  =  N  +•  1  ; 


JU, 

t:ii 


i 


BEGIN 


£ 


V- 


r* 

,v 


NewDa ta  C  K3.  Org  i n  -  St  I  d  ent  C  J3.  Org  i  n ; 

K  :  *  K  -  1 ; 

END; 

N  :  =  1 ; 

K  :  =  S  taNu  mb ; 

WR  ITELNdrai  LOut) , 

UR  ITELNdrai  lOut)  ; 

WR  ITELNdrai  lOut), 

WR  ITELN  dr  ai  lOut/  dhe  two  stations  par  t  i  c  i  pat  i  ng  in  the  intersection  "), 
WR  ITELN(Trai  lOut,  'calculations  are:  ’  ) , 

WR  I TELN  drailOut,  '  Station#',  Ne wData CN3  .  Org  in:  1  ) , 

WR  ITELNdrai  lOut,  '  Station  # ',  NewDataCN  +  13.  Orgin:  1), 

WR  I  TELN  drailOut)/ 

WR  ITELNdr  ai  lOut,  dhe  following  information  determines  the  back'), 

WR  ITELNdrai  I  Out,  '  LDB  "  s  to:  ')■ 


A  :=  NewDa  taCNJ.  Or  g  in; 

13  ;=  NewDa t a IN  +  13.  Org  in; 

A  :  =  NewDa ta C A 3 . Or g  in; 

B  :-  NewDa  taCB  3.  Org  in; 

CASE  S  taNu  mb  OF 


2  :  WRITELNC  Trai lOut,  ' 

3  WRITELNv Trai lOut,  ' 


No  other  station'); 

St  a  t  ion  # ' ,  NewDa  ta  C  K  3 .  Orgin:  l)i 


4  :  BEGIN 

WR  I  TELN  (Trai  lOut,  ' 
WR  ITELNdrai  lOut,  ' 

END; 


Station  #  ' «  NewDa  t  a  CK3 .  Or  gin:  1 
Station  #',NewDataCK  —  13.  Org  in 


l\* 

v> 


5  :  BEGIN 

WR  ITELNdrai  lOut,  ' 
WR  ITELN(Tra  i  lOut,  ' 
WR  I  TELN  drailOut,  ' 

END, 

END, 


Station  #  ' ,  NewDa  t  a  CK  3 .  Orgin:  1  )  . 
Station  #  ' ,  Ne  wDa  t  a  CK  —  13  Or  gin 
Station  #',NewDataCK  —  23.  Orgin 


WR I  TEL N(TrailOut) , 

WR  I  TELN  drailOut,  ' - 

WR  I  TELN  (  Tr  a  l  1  Du  t ,  ' - 

WR  ITELNdrai  lOut,  ' - 

WR  I TE ( Tr a  i  lOu t ,  ' - 

WR ITE ( Tra l  lOut ,  New Data  CN 
WR  I  TELN  (  Trail  Gut), 


- '  )  , 

- '  )  ; 

- Beginning  Of  Inputted  Data - 

-For  Intersecting  Station" s  ' , NewDat a  CN J .  Or g i n :  1 , 
-  1  3 .  Org  i  n :  1 ,  ' - '  )  ; 


WR I TELN ( Tr a i lOu t ,  ' - )  > 

WR  ITELN(Trai  lOut.  ' - ')) 

WR I TEL  N(TrailOut) ; 

WR  I  TEL  N(TrailOut); 

WR ITE ( Tra 1 1  Qu  t ,  'Intersecting  coordinates  for  Station  #  NewDataCNl .  Qrg  in.  1  ’ 
WR  ITE(  Trai  lOu  t ,  '  and  Station  #  NewDataCN  +  13.Qrgin:l»'  is:'); 

WR ITELN (TrailOut) , 

WR ITELN (TrailOut) ; 

WR  ITELN  (TrailOut,  '  (  ' ,  In  t  x  .  6 :  3,  ' ,  ' .  Inty :  6:3, 

WR ITELN (TrailOut) ; 

SupportCount  :=  0; 

FOR  J  =1  TO  StaNumb  DO 

IF  < (StI dentCJl.  Ref  lob  =  TRUE)  AND  < St I d ent C J3. F  lag  =  FALSE))  THEN 
BEGIN 

Allstate  J3.  Flag  :=  S 1 1  d  e  ntC  J  3  .  F 1  a  g ; 

WR ITELN < TrailOut), 

WRITELNI Trai lOut,  'X-c oor d ina t«  ' ,  St  I d en t C J1 .  Xcoo r d : 6 : 3,  '  of  Station  F 
StI  dent  CJ1.  Qrg  in:  1)  i 

WRITELM( Trai lOut,  'loaded  into  the  system.  '); 

WRITELN( Trai lOut ) ; 

AllStatLJl.  Xcoard  :  =  StldentCJI.  Xcoord; 

WRI.TELN<  Tra  i  lOut,  'Y-coor  d  inate  ' ,  St  I  d  en  t  C  J1 .  Yco  o  r  d  :  6  :  3,  '  of  Stacien  -*■ 
St  I  den t C J1 .  Org i n :  1 ) ; 

WR ITELN  <  Trai lOu  t ,  'loaded  into  the  system.  '); 

WRITELN( TrailOut); 

AllStatt  Jl.  Vcoord  :  *=  St  I  d  en  t  CU3 .  Vc  oo  r  d ; 

DistX  :=  (Intx  -  A1 1 S ta t C J3 .  Xcoo r d > ;  (*  could  be  ABS  *) 

DistY  : =  (Inty  -  A1 1 Stat C J1 .  Ycoord ) ;  (*  could  be  ABS  *) 

IF  DistY  =  00.000  THEN 

BEGIN 

WR I TELN ( TrailOut ) ; 

WR ITELN ( Trai lOut , '  This  station  can  not  exist  at  the  location'): 
WRITELN(  Trai  lOut ,  '  of  (  DistX:  6:  3.  ',',DistY:6:3.  ')'); 

WR11 ELN( Trai lUut  ) ; 

END 

ELSE 

BEGIN 

IF  ((DistY  <  00.000)  AND  (DistX  <00.000))  THEN 
BEGIN 

AngleTheta  :  =  AR CTAN ( Di s t X/D i stY )  *  180/Pie, 

AngleTheta  :=  AngleTheta  +  1Q0.  QQQ; 

IF  AngleTheta  <  00.000  THEN 

AngleTheta  :=  360.000  +  AngleTheta; 

END 

ELSE  D-1.4 


IF  < (DistY  <  00.000)  AND  (DistX  >  00.000))  THEN 
BEGIN 

DistY  :  =  00.  000  -  DistY, 

AngleTheta  :=  ARCT AN< D i st X /Di s tY )  *  180/Pie; 

END 

ELSE 

BEGIN 

AngleTheta  : =  ARCT AN( Di stX /DistY)  *  100/Fie, 

END; 

WRITEUM( Trai lOut ) , 

WRITELN< Trai lOut,  'Th e  exact  computed  LOB  angle  is  =  ' , Ang leTheta .  r 
WRITELN<  Trai lOut  ) , 

Temp4  :=  ABS (StI dent CJD . Theta  -  AngleTheta), 

Stld ent C JD. MissAngle  :  —  Temp4; 

Allstate JD.  MissAngle  :=  StldentCJD.  MissAngle; 

IF  <  T emp  4  >=  360  O)  THEN 

Temp4  :=  Temp4  —  360.0) 

IF  ( T  emp  4  <=  (3  *  Sigma))  THEN 

BEGIN 

Temp5  :  =  Temp4  *  TempA; 

Tab  1  eCA,  BD.  Support  :=  TempS  +  Tab leC A, BD . Sup p ort ; 

Tab  1  eCB,  AD.  Support  :  =  TempS  +  Ta b le C 8. AD . Sup p or t ; 

Tab  1  eCA,  BD.  SupportFlagCStldentC  JD.  Org  inD  :=  TRUE; 

TableCB,  AD.  SupportFlagCStldentCJD.  OrainD  :  =  TRUE; 

SupportCount  :=  SupportCount  +  1; 

Tab  1  eCA,  BD.  SupCountFIag  :=  SupportCount; 

END 

ELSE 

BEGIN 

Tab  1  e  C A,  B  D.  S upp  o r  tF  1  ag CS 1 1 d  e n t C  J D .  Or g  in  D  :=  FALSE; 

TableCB,  AD.  SupportFlagCStldentCJD.  Or  g  i  n  D  ;  =  F  ALSE ; 
WRITECTrailOut,  'Station  #  St  Id  entC  J  D .  Or  g  in  :  1 ,  '  "  s  exact  LCB  '  )  , 
WRITE(Trai  lOut,  '  is  greater  than  3  sigma. 

WRITELN< Trai lOut ) , 

END, 


WR ITELN ( Trai lOut  ), 

WF I TELN ( Tra l lOu t , 'Difference  between  the  back  LOB  angle  and  the'), 
WRITE  (Trailout,  'actual  angle  for  Station  #',StIdentCJD.Qrgin:l>, 

WR  ITEdrailOut,  '  =  '  ,  Temp  4  6  :  3,  '  .  ' ) ; 

WR ITELN ( Trai lOut  )  ; 

WRITELN  (  Tra  i  lOut  )  ,  D--15 


J* 

ID 

u 


;■  l 


jr 


i 


END; 


END 

ELSE 

IF  <  (StldentC  J3.  Ref  lob  =*  TRUE)  AND  (StIdentCJ3.  Flag  »  TRUE))  THEN 
BEGIN 

AllStatCJ3.  Xcoord  :  -  StIdentCJ3. Xcoord; 

Al lStatCJ3 . Ycoord  :«  StIdentCJ3.  Ycoord; 

Al lStatrJ3.  Alpha  StIdentCJ3.  Alpha; 

Al lStatCJ3. Theta  :=  StIdcntCJ3.  Theta; 

Al lStatTJ]. Ref  lob  :=  St  I  den t C J3 .  Ref  lob; 

Al lStatCJ3.  Flag  :  =  St IdentCJ3.  Flag; 

Al ISta t C J3 . Mi s sAng 1 e  :m  St  I  den t C J 3 .  Mi s sAng 1 e; 

A1  lStatCJ3.  Org  in  :»  St  I  den tC J3 .  Org in; 

END 

ELSE 

BEGIN 

Al lStatCU3.  Xcoord  :  =  00.  OOO; 

Al lStatCJ3. Ycoord  :  =  00. 000; 

Al lStatCJ3.  Alpha  :  -  00. 000; 

A1 lStatCJ3.  Theta  : *  00. 000; 

Al lStatCJ3. Ref  lob  :»  FALSE; 

Al  lStatCJ3.  Flag  :»  FALSE; 

Al  IStatCJl.  MissAng  le  :■  00.000; 

Al IStatC J3 . Org in  :■  St  Iden t CJ3 .  Or g  in; 

END; 


WR  I  TEL  N(TrailOuti  ' - 

WR  ITELN(TrailGut,  ' - End  Of  Inputted  Data - '), 


WR ITE( Trai lOut,  ' - For  Intersecting  Station"s  ' > NeuData CN3 .  Org i n :  1 , 

WR  ITE  (  Tra  i  lOut ,  NeuiData  CN  13.  Orgin:l.  ' - '); 

WR ITELN( TrailOut) ; 

WR  ITELN(Trai  lOut.  ' - ' )  ; 

WR I TELN (TrailOut); 

WR ITELN(TrailQut) ; 

WR ITELN( TrailOut) ; 

WR  I  TELN  (TrailOut); 

WRITELN(TrailOut); 

WR ITELNCTr ailQut) ; 

WR  I  TELN  (  '  Leaving  ILOBS'); 

TempD  ;  =  00.  000; 

END;  (*  PROCEDURE  "ILOBS"  *) 

(  ++#**#++**#  +  ##+#***  *******#******#**###*#■*■»•**■**■*•#•**■»•■#■*■*■»••»•*••«■  * 

PROCEDURE  XSEC  (  Stldent  :  Lobtype; 

StaNumb  .  INTEGER; 

I  :  INTEGER; 

J  INTEGER; 

VAR  I  ntx  REAL; 

VAR  Inty  :  REAL;  0-16 

VAR  Computed  BOOLEAN); 


I  (****^***** ***-»***#*•»-»**•►***••♦**•****♦♦■»*■****♦***•»♦* *********HHHi'*****'iH*-** *•**•** 
(* 

This  procedure  uses  the  two  LOBS  to  determine  if  there  is  an 
intsrsecton  using  the  criteria  for  a  valid  fix  estimate, 
intersecting  coordinates  are  INTX  and  1NTY.  "  XSEC“  establishes  a 
series  of  LOBS  called  INTLO  from  all  stations  to  the  intersection 
point. 

*) 

VAR 

■y  X,  Y 
H  X  d  i  st 
Yd  i  st 

^  Listl,  List2 
Y  AllStat 


BEGIN 


WRITELN('  Entering  XSEC'); 


X  :  =  I ntx; 

Y  : =  Inty; 

Xd  ist  :  *  00.  000; 

Yd  ist  :  »  00.  OOOi 

Intx  :  *  00.  OOOi 

Inty  :  =  00.  000; 


Xdist  : =  (StldentC 13.  Xcoord  -  Stld entC J3.  Xcoord ) ;  (*  could  be  ABS  *5 
Ydist  :»  (StldentC  13.  Ycoord  -  StldentC  J3.  Ycoord  >;  (*  could  be  ABS  *) 

XCOhff5  (  I  >  Ji  Xdist.Ydist.  Stld  ent.  Int  x  »  Inty.Computed); 

IF  (Computed  =  FALSE)  THEN 

BEGIN 

WR I TELN (TrailOut); 

WR I TELN  (TrailOut.  '  No  intersection  found  from  this  dat.--.,  using'  ), 

WR  ITE (  Tra  i  lOut  •  '  Station  #'.StIdentCI3.Qrgin:l/'  and  '  ) ; 

WR ITE ( Tra i lOut «  '  Station  # ' . St IdentC JI . OTg in: 1 > ; 

WR I TELN (TrailOut) ; 

WRITELN(TrailOut); 

END; 

i  WR ITELN ( '  Leaving  XSEC'); 

. END-  <*  PROCEDURE  "XSEC"  *) 

'  «-*  *  *«-*#  **■»*•»  <H****W***#**^-*#***  •#.■**■*  **■«•■*  «««*##«  \ 

PROCEDURE  ZEXMVAR  Table  ;  TVT ; 

|  Lob  Lob  type) ;  d-17 

(*  This  procedure  checks  the  validity  of  the  intersections.  If  it  fails 


REAL;  <*  Dummy  variables  for  Intx  and  Inty  *> 

REAL;  (*  Distance  of  XI  -  X2  *) 

REAL;  (*  Distance  of  Y1  -  Y2  *) 

Intlotype;  (*  used  for  computing  intercept  coords  * ? 
Lobtype;  ( *  data  structure  that  will  contain  all  of 
the  LOB's  that  return  to  station  that 
did  not  attempt  any  LOB's.  *) 


the  test  for  validity  <  i .  «.  >  0.  b  )  than  tha  corresponding  value  in 

tha  Truth  table  will  be  marked  indicating  that  tha  intersection  failed 
the  validity  test. 

*) 

VAR 

BinNumber  :  INTEGER; 

AStation  INTEGER; 

BEGIN 

WRITELNI'  Entering  ZEXI'); 

FOR  BinNumber  :=  1  TO  5  DO 
FOR  AStation  : =  1  TO  5  DO 
IF  Tab  1 aCBi nNumber » ASta t ion  3.  FI ag  =  TRUE  THEN 

IF  <ABS(LobCBinNumber  3.  Xcoord  -  LobCB  inNumber  3.  Yc  oor  d  )  )  >0.6  THEN 
Tab l e CBinNumb er»  AStat  ion3.  Flag  :■  FALSE; 

WR  ITELN<  '  Leaving  ZEXI'); 

y 

END;  <*  PROCEDURE  "ZEXI"  *) 

( *•»*■»*■»•**■»***■»**•»•■*#*****•»*#*•»■***•****#***•**#***■**■*****'«-»#**((-»•»•»**■«■■«■»■»■»■»■*•»*•»■**-<• 
PROCEDURE  RECONVERT (  Lobrec  :  Lob  type; 

VAR  Starec  :  Lob  type); 

(  ■»•**  ■*■***  •*#•**  #**■#  -»**-»  *#*«  **•»*  *•»*■»  *■**# -»******# ************  ****-»*-»«  #****"***'!" 

(■*  This  subroutine  takes  the  data  stored  in  "Lobrec "and  places  it  in  the 
original  order  in  the  data  structure  called  "Starec". 

*) 

VAR 

J  :  INTEGER; 

Location  :  INTEGER; 

BEGIN 

WR ITELN<  '  Entering  RECONVERT '); 

FOR  Location  :  =  1  TO  5  DO 
BEGIN 

U  :■  LobrecCLocation].  Org  in; 

StarecCJ].  Xcoord  :  ■  LobrecCLocation].  Xcoord; 

Star  ecC  J].  Yc  oord  :*  LobrecCLocation].  Ycoord; 

Star ecCJ3.  Th eta  :®  LobrecCLocation].  Theta; 

StarecCJ3.  A1  pha  :  =»  Lobr  ecCLocat  ion].  Alp  ha; 

Star  ecC  J3.  RefLob  :»  Lobr  ecCLocat  i  on  3.  Ref  Lob; 

Star  ecCJ3.  Flag  :■  LobrecCLocation],  Flag; 

END; 

WR  I TELN  (  '  Leaving  RECONVERT'); 

END;  (*  PROCEDURE  "RECONVERT"  *) 

0.18 

<  ***********************************#■****#**###**#***#*#******#**.****#**  *.#*.*.(*.* 
PROCEDURE  CONVER T (  Starec  :  Lobtype; 


£ 


i 


VAR  Lobrec  :  Lobtype); 

(*  This  subroutine  takes  the  data  stored  in  "Starec"  and  places  it  in  a 
sequenced  order  starting  in  the  first  and  subsequent  cells  of  the  data 
structure  "Lobtype". 


VAR 


e 


J,  K 

Locati on 


BEGIN 

WRITELNi' 


INTEGER ; 
INTEGER; 


Entering  CONVERT'); 


e 

O. 


f<< 

& 


J  =  0. 

K  =6, 

FOR  Location  : =  1  TO  5  DO 

IF  S tare c CLo cat i onl.  Ref  1 ob  *  TRUE  THEN 


i 


$■ 


BEGIN 

J  :  =  J  +  1;  ( *■  i  nde 

Lubr  et-t  J  ].  Aboard 
Lobr  ecC  J3.  Yc  oord  :  = 
LobrecC  J3.  Th  eta  ;  = 
LobrecCU3  Alpha  ;= 
LobrecC  JI.  Ref  lob  ;  = 
Lobr ec C J 3 . FI ag  .= 
LobrecC J3.  Orgin  ;  = 


for  Lobrec  *) 

S  Ca  recCLocat ion 3 . Xcoord; 
StarecCLocationl. Ycoord; 
Star  ecCLocat  ion].  Theta; 
Star ecCLocat ionl . Alpha; 
StarecCLocat ionD. Ref  lob; 
Star ec CLocat i on3 .  Flag; 
Loca  t ion ; 


END 

ELSE 


BEGIN 


5 

iV 

M  * 

V- 


■„*> 


LohrecCK3.  Xc  oord  ;  =* 
LobrecCKI  Ycoord  :  =* 
Lobr ec CK 3 . Th e ta  : — 
LobrecLK],  Alpha  ;  = 
Lobr ec CK 3 . Re f  lob  ;  = 
Lobr ec IK 3 .  FI ag  :  = 
Lobr  ec  CK I .  Or  g  in  :  =* 


for  Lobrec  *) 

Star ec CLocat i on 3 .  Xcoord; 
Star ec CLocat i on 3  .  Ycoord; 
Star ec CLocat i on 3  Theta; 
Star  ecCl.ocat  ion3  .  Alpha; 
Star#cCLocation3. Reflob; 
Star ecCLocat i on  3 . Flag; 
Loca  t  ion ; 


END; 


s.  WR I TELN  <  '  Leaving  CONVERT'); 

til 

END,  (+  PROCEDURE  "CONVERT"  *> 


c  ** ■******♦*■» a*************** * *** a*******  a*#* *■»■*******«-»* a-*** *****<»**4»* * *■ 
PROCEDURE  GAXKVAR  Table  ;  TVT ; 

VAR  Computed  BOOLEAN; 
i*  Dspxi  :  INTEGER; 

“  Starec  Lobtype; 

_ StaNumb _ :  INTEGER), 


D- 19 


<  **i1- ****  ***#»■*•**  *■*■«■#•  *#*■«■»*■**•*-»-*■*  ■*■*** *#*-)M»'*** 

<*  This  procedure  extracts  the  data  from  the  database  for  use  in  iC 

computing  the  intersection  of  the  LOB  '  s  from  the  DINSETS  for  each 
station.  Parameter  "Table"  will  contain  truth-value  assignments  for  u\, 

valid  intersections  in  the  Einsets  of  each  Station. 

* ) 

VAR 


1/  J 

INTEGER 

B i nNumb  er 

INTEGER 

AStat i on 

INTEGER 

Ne  x  tCri  c 

INTEGER 

Ma  xSuo  port 

INTEGER 

Temp Station 

INTEGER 

La opCn  tr 

INTEGER 

Locationl 

INTEGER 

Vo  cat l on2 

INTEGER 

Lo  c  at i on3 

INTEGER 

Lo  cat i on4 

INTEGER 

XI 

REAL; 

Y 1 

REAL; 

inc* 

REAL; 

Inty 

REAL; 

MinExactLob 

REAL; 

Lo  b  r ec 

Lob  type, 

All Sta  t ion  s 

Lob  type 

(*  index  for  the  data  structure  LOBREC 


(*  keeps  track  of  the  back  LOB's  to  the 

stations  whose  intersection  has  not-  yet 
been  calculated.  *) 


BEGIN 

WR  I TELN  (  ' 


Entering  GAXI  ' ) ; 


AStat ion  : -  2 ; 
TempStation  :=  AStation; 
B i nNumb  er  : =  1 ; 

In  tx  ;  =  00  000; 

Inty  ;  =  00. 000, 


CASE  StaNumb  OF 


END; 


Loop  Cntr 

:  =  1; 

Loop  Cntr 

:=  3; 

<k„  » 

Loop  Cntr 

:  =  6, 

Loop  Cntr 

:  =  10; 

c- 

T  * 

V' 

T (Starec , 

Lobrec  ) ; 

E? 

Dsp  x i  <= 

LoopCntr  DO 

BEGIN 

IF  Temp  Station  =  StaNumb  +  1  THEN 
BEGIN 


AStation  :=  AStation  +  1; 
Temp  S  ta  1 1  on  :  =  AStation; 

B  l  n  N  u  m 


i 

l  ■  . 

k 

ux 

tt 

if* 

i  •. 

L'- 

t 


& 

it- 


K 


Lob r e c C B inNumber  -  13. Flag  :  *  FALSE; 

END; 

IF  ( ( Lob r ec C B inNumb er 3 . Re f 1 ob  *  TRUE)  AND 

( Lob r ecCTempStat i on3  .  Ref 1 ob  ■  TRUE))  THEN 

BEGIN 

XSEC  <  Lob r  ec,  S taNumb ,  B  i  nNumber .  TempStation,  Intx,  Intiji  Computed)  ; 
END; 

Locations  :  *=  Lob r e c CB i nNumber 3 .  Or g in ; 

Location4  :=  Lab r e c C T omp S tat i on  3 .  Or g i n; 

IF  Computed  =  FALSE  THEN 

BEGIN 

TableCLocation3,  Locat  ion43.  Flag  :  *  FALSE; 

Tab  1  e  CLoc  atian4,  Loca  t  ion3  3.  Flag  :  =  FALSE; 

END 

ELSE 

BEGIN 

Tab  1  e  CLoc  at  ion3,  Locat  iar.4  3 .  X  inter  cep  t  ;=  Intx; 

Tab  i e  CLoc ation3.  Locat  ion4 3.  Yi nter  cep  t  :=  Inty; 

Tab  1  e  CLoc  at  ion4,  Locat  ion3  3.  X  inter  cep  t  :  *  Intx; 

Tab  1 e CLoc at i a n4,  Locat ion3 3 .  Yi nter cep t  .=  Inty; 

Tab  1  e  CLoc  a  t  ion3,  Loca  t  i  on4  3 .  FI  ag  :=  TRUE; 

Tab  1 e CLoc at i on4. Locat i on3 3 . FI ag  :*  TRUE; 

LcbrecCBinNumber 3 . Flag  ;  *  TRUE; 

Lobr e c CTempS t a t i on 3  Flag  : ■  TRUE; 

ILOBS < A1 IStat ions > Lob rec <  Intx,  In t y , St aNumb , Tab le ) ; 

FOR  I  :=  1  TO  StaNumb  DO 

IF  AllStati  onsC  13.  Flag  =  FALSE  THEN 
BEGIN 

J  :  =  LobrecCIl.  Orgin; 

Star  ec  C  03.  M  issAng  1  e  :=  A1 IS tat i ons C  I  3 .  M  i  ssAng 1 e , 

END. 

END; 

Dspxi  :  =  Dspxi  +  li 
TempStation  : —  TempStation  +  1, 

Lobr  e  c  C  TempSt  a  1 1  q  n  -  13.  Flag  :=  FALSE; 

END, 

RECONVERT  (Lobrec.,  Starec), 
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WR I tElN ( TrailOut) . 


WRITELN  (TrailOut,  'This  it  th*  table  that  shows  th*  relation  of  data')/ 

WR  ITELN(Trai  lOut,  'between  ang  of  th*  station*.  ' ) ;  fc  { 

WRITELN(TrailOut) ; 

WRITE <  Trai  lOut,  ' -  i 

WRITE(Trai lOut.  '  ');  1 

WR ITELN(Trai lOut ) ;  9 

WRITEtTrai lOut.  'Int*r**cting  Station*  Xintercept  Vinter 

WRITEITrai lOut,  '  Support');  > 

WR ITELN(TrailOut) *  ( 

WRITE  (Trai  lOut,  ' -  -  - 

WRITE  (Trai  lOut,  '  '); 

WRITELN(TrailOut) ; 

MinExac tLob  lOOOOOOOOOO. 000 ; 

Ne  x  tOn  e  :  »  5; 

MaxSupport  :=  0, 

FOR  BinNumber  :  =  1  TO  5  DO 
FOR  Temp Station  : *  1  TO  5  DC 

BEGIN 

WRITE(TrailOut,  'Station  BinNumbtr:  L  '  Station  #'»  Temp  Stat  i  on :  i  ) ; 
WRITE(TrailOut>  '  '); 

WRITE  (Trai  lOut,  Table  CBinNumb  *r»  T  empStat  i  on3.  X  intercept:  6:3); 

WRITE  (TrailOut,  '  '); 

WRITE(TrailOut.  Table CBinNumb er*  TempStat  i  on3.  Yintercept:  6:  3); 
WRITE(TrailOut,  '  Tab loCBinNumber,  TempStationl.  SupCour. tF  lag  6  : 

WRITELN( TrailOut); 

IF  ( ( (Tab leCBinNumber, TempStation3.  SupCountFlag  <  NextOne)  AND 

(Tab leCBinNumber, T*mpStation3.  SupCountFlag  >  Ma xSup p or t ) )  AND 
(Tab leC BinNumber*  TempStati on].  Flag  -  TRUE))  THEN 

BEGIN 

MaxSupport  :*  Tab leCBinNumber*  TempStati  on3.  SupCountFlag; 

Locationl  :  =  BinNumber; 

Location2  :»  TempStation; 

END; 

END; 


FOR  I  : -  1  TO  5  DO 
BEGIN 

WR  I TELN  (TrailOut); 

WR  I TELN (TrailOut*  ' _ '); 

WR  ITELN(TrailOut,  'Station  1:1,  '  data:'); 

WR  I  TELN  (TrailOut.  ' - '  ); 

WR  I  TELN  (TrailOut); 

WR  ITELN  <  Tr  a  i  lOu  t«  'X-coord  -  ',Star*cCI3.  Xcoord.6:3); 

WR  ITELN<Tr  ailOut*  'Y-coord  «  '*  Starec  C  13.  Ycoord:  6:  3) ; 

WRITELN(TrailOut,  'Theta  »  ',  Starec  C  13.  Thet/i:  6:  3); 

WR  ITELN(  Tr  a  i  lOut,  '  Alph  a  *  Starec  CI3.  Alpha:  6:  3); 

WR  ITELN'Trai  lOut,  'Ref  1  oh  -  ' , Starec C  1 3 ,  Ref  1 ob ) ; 

WR  ITELNdrailOut,  'Flag  D-22  »  ',  Starec  C  13.  Flag  ) ; 

WR ITELN( Tr a l lOut,  'Semimajor  -  ' , Starec C  1 3 ,  Semi ma jor : 6:  3 > ; 
WR  ITELN(Trai  lOut,  'Semiminor  -  '»  Starec  C 13.  Semi  minor :  6:  3 ) ; 


WR  ITELN  (TrailQut.  'Distance  ™  '.StarecEI].  Distance:  6:  3); 

WR  ITELN ( Tr  a  i lOut.  'Qrgin  =  ' » Starec C 13 .  Org i n: 2) ; 

WR  I TELN (TrailOut). 

END. 

WR  ITELN(Trai  IGut,  'From  the  data  submitted  by  each  station  and  the  data  ). 

WR  ITELN(Trai  IDut,  '  extracted  from  various  r.alcu  lati  ons  the  best  "FIX"'); 

WR ITELN(TrailOut,  '  location  of  the  object  attempting  to  transmit  is  at  '); 
WRITELNdrailOut); 

WR ITE< Trai lOut.  '  (  ' ,  Tab letLoca t ion  1 , Lo cat i on23 . X  in t er c ep t : 6: 3.  '■ 

WR  1TE<  Trail  Out.  TableCLocationl.  Loc  at  i  on21 .  Yintercept:6:3»  '  )  ')i 
WR  I  TEL  N(TrailOut); 

WR ITELN(TrailOut) ; 

WR  I  TEL  N(TrailOut), 

WR I  TEL  N(TrailOut).  ' 

WR ITELN (TrailOut); 

WR I TELN (  '  Leaving  GAX I  ' ) ; 

END.  <*  PROCEDURE  "GAXI"  *) 

(  *-»■*****  **•**■*-*•**••*-*■*-*  ***-»■*-**•*-«■*-*■*  -**-**•****■*■■*■«*  *■»**  ■*■»**■»■-»** -a-**-**^- 

PROCEDURE  SUXI (VAR  Dspxi  :  INTEGER); 

(  *-**  ««-»* *****-»**•»-»** **4**«a***xa**  **-**■■**■*•«■■*■**■«■  +*#+**-***+*■*#*-***#*■# ■*■*■*■»  »•*•*■*  -is 

(4  This  procedure  sets  up  the  d i p lac em.en t  variable  (  parameter  Dspxi  ) 
and  sets  up  the  ability  to  extract  data  from  the  database  through 
use  of  the  displacement  variable. 

*  • 

BEGIN 

WR I TELN ( '  Entering  SUXI'); 

Dspxi  =  1;  ( -*  -1  in  actual  program  *) 

WR I TELN ( '  Leaving  SUXI'); 

END,  <*  PROCEDURE  "SUXI"  *) 

4*4*  *4  4  **W  *  ****  4 -*-*■»■  *■*■■»■*  **-**•  ■*■#■-»•  **-4-*-*  *******  *************  *******************  4*4 

PROCEDURE  LMARtX  ( VAR  Table  .  TVT ; 

Lobrec  Lobtype  ); 

(  ***tHHS-#**#***#<HHHt*m**#**»4**«4***4*H**-*»*****#*#**t*1Ht**««4<HHHt'iHtiHH-»i;-tt-H 

(4  This  procedure  marks  LOBS  acceptable  for  use  in  fix  estimation  Marked 

LOBS  for  optimixed  fix  estimation  are  used  by  "FPEAK"  to  obtain  a  best  fix 

*) 

VAR 

knnNumoer  INTEGER; 

I  :  INTEGER, 

DEGIN 

WR I TELN  <  '  Entering  LMARK ' ) ; 

I  =  1  , 


FOR  QinNumber  . =  1  TO  5  DO 
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»?i  Wfi *w»  #1  in  p  r\  nr  ihiwpj i^mihi 


IF  ( ( Lobr ecCE inNumber 3.  Xc oord  *  Lobrec CD inNumber  +  ll.Xcoord)  AND 
(LobrecLB inNumber].  Ycoord  =  Lobr ecCB inNumber  +  13.Ycoord)>  THEN 
Tab l e CAStati on,  B inNumber 3. Flag  :  =  TRUE; 

I  : **  I  +  1  ; 


WR  ITELN< 


Leaving  LMARK'); 


END;  <*  PROCEDURE  "LMARK"  *> 


PROCEDLIRE  FI  RSTP  ASS  ( Lobr  ec  :  Lob  type;  ± 

StaNumb  :  INTEGER;  m 

Table  :  TVT ;  ^ 

Intlo  :  In  t la t  y  pe ) ; 

(  #****^**^*HHHt******^**^**********^**#**#****-*-H-**-*-»*-S-3--ft**'*tf  *  *-*#*-**'*0 

(  *  This  sub  does  whatever. 


Dspxi  :  INTEGER; 
Computed  BOOLEAN; 


BEGIN 


WR ITELN< 


Entering  FIRSTPASS'); 


Computed  : =  FALSE; 

SUXI (Dspxi); 

GAXI<Table, Computed, Dspxi,  Lobrec,  StaNumb ) ; 


WR  ITELN<  ' 


Leaving  FIRSTPASS'); 


END,  (*  PROCEDURE  "FIRSTPASS"  *) 

PROCEDURE  LOADDATA( VAR  Lobrec  ;  Lobtype; 

VAR  StaNumb  ;  INTEGER); 

<*  This  procedure  prompts  the  user  for  input  to  the  TRAILBLAZER  program. 
*> 


Te  mp6 


BEGIN 


INTEGER; 


StationCntr  :  INTEGER; 


WR ITELN( ' 


Enterina  LOADDATA'); 


REPEAT 


WR I TELN ( Tra 1 IDut ); 

WR I TELN < Tra i lOut ,  '  How  many  stations  will  be  reporting  LGB"s') 
WR ITELN ( Tra i lOut ,  '  on  the  transmitting  object.  '); 

READ (Tra i 1 In, StaNumb ); 

WRITELN(TrailOut, StaNumb );  D-24 


Temp  6  :■  StaNumb; 

IF  <  (Tempi*  <«  0)  OR  (Tempi*  >«*  6)  )  THEN 
BEGIN 

WR1TELN(  Trai lOut ) ; 

WRITELN< 'This  value  ' ,  StaNumb :  2>  '  is  invalid'); 

WRITELN( 'Try  again  with  a  valua  from  1  to  3.  '); 

WRITELN(  Trai  l Out ) ; 

END; 

UNTIL  ((Temp6  >  0)  AND  (Tempi*  <  6) >; 

REPEAT 

BEGIN 

Tempi*  :  =*  Tempi*  -  1; 

REPEAT 

REPEAT 

WRITELN(TrailOut); 

WR  ITELN(TrailOut»  'Input  tha  Station  that  is  "FIXING".  '); 

READ(Trai 1 Inj  StationCntr); 

WR  ITELN(TrailOut»  StationCntr); 

IF  ((StationCntr  <«  0>  OR  (StationCntr  >«  6))  THEN 
BEGIN 

WR ITELN(TrailOut)  i 

WR  ITELN(TrailOut»  'This  valua  '•  StationCntr  :  2*  '  is  invalid'); 
WRITELN(TrailOut«  'Try  again  with  value  from  1  to  5. 

WR ITELN(TrailOut) ; 

END; 

UNTIL  ((StationCntr  >  0)  AND  (StationCntr  <  6)); 

IF  Lobrec CStationCntr 3.  Ref lob  *  TRUE  THEN 
BEGIN 

WRITELN(TrailOut); 

WRITELN(Trai lOut,  '  Station  StationCntr: 1>  '  has  already  been'); 
WRITELN( Trai lOut, '  processed.  Try  again  with  a  new  station.  '); 

WR I TELN ( Tra i 1 Ou  t ) ; 

END; 

UNTIL  LobrecCStationCntr3.  Ref  lob  ■  FALSE; 

WR  ITELN(TrailOut) ; 

WR ITELN(Tr  ailOut*  '  Input  the  x-coordinate  of  Station  # ' » Sta t ionCntr :  1 ) ; 
RE  AD  (Trail  In,  Lobrec  CStationCntr  3.  X  coord  ); 

WR  1TELN (Tr ailOut*  Lobrec CSt at i onCnt r 3.  X c oor  d :  6:  3); 

WR ITELN(TrailOut) ;  D-25 

WR !TELN(TrailOut,  ' Input  the  y-coordinate  of  Station  StationCntr : 1 > ; 


READ<TrailIn,  Lobrec  CSta  t ionCntr  ].  Vcoord); 

WR  ITELN  (TrailOut,  Lobrec  CStati onCnt r  3.  Yc  oor d : 6:  3); 

WR  I TELN (TrailOut); 

WR  ITELNi'Trai  lOut,  'Input  the  angle,  in  relation  to  true  North,  formed' 
WR  I  TELN  (TrailOut.  'by  Station  #  ' .  St  a  t  i  onCn  t  r :  1  <  '  and  the  possible  '); 
WR  I  TELN (TrailOut,  'location  (  FIX  )  of  the  transmitting  object'); 

READ< T rail  In,  LobrecCStationCntrl.  Theta); 

WR  ITELN  CTr  ai  lOu  t,  LobrecCStationCntrl.  Theta:  6:  3); 

WR  I  TELN  (TrailOut); 

IF  Lob r ec CSta t i onC ntr 3 .  Theta  >  ISO.  00  THEN 

LobrecCStationCntr3 . Theta  :  =*  Lo b r ee  CSta  t  ionCntr  3 .  Th  e ta  -  360.00; 
Lobrec CStationCntr 3. Al pha  90  -  LobrecCStati onCntr3.  Theta; 

Lo brec CSta t ioncntr 3 . Re f lob  :=  TRUE; 

WR  ITELN  (TrailOut); 


UNTIL  (Temp6  =  0), 

WR  ITELN (  '  Leaving  LOADDATA'); 

END;  <*  PROCEDURE  "LOADDATA"  *) 


(  ##*-tt-#**#*****-»-a-*#**-#^*-fr  ■«■■«■■**  #■#**■**#■«■  **■*•#■#■*■*■#■■**■**  *«"»-»****«#-**  •«•**#«•■#■*■*  -fHHH 


PROCEDURE  ESTMP (Lobrec 

IntPass 
Index 
I  nt  1  o 
Table 


Lobtype; 

INTEGER; 

INTEGER; 

Intlotype; 

TVT); 


###•#  «  ■*  > 
<V 


(#*#  ***■&  #■»**•»****«**  ***■■*  ********  #**••*■**■*•«•  ****  ********  ********  ************  ** 

<* 

This  procedure  is  called  by  "HPFIX"  to  obtain  a  best  fix  t 

estimate.  A  single  fix  or  one  of  mutilple  fixes  in  either 
the  normal  mode  or  the  degraded  mode  may  be  requested. 

"ESTMP"  calls  the  following  procedures  and/or  functions: 


SUBS 


COMMENTS 


1 ) 

GAX  I 

Ex  t 

racts  tw 

O 

r 

nes 

of  Berin 

com 

puting  i 

nter 

sed 

t  i  on. 

2) 

INELL 

To 

determin 

e  an 

d  ensure  tha 

doe 

sn't  dup 

1  ica 

te  the  last 

3) 

LMARK 

Mar 

ks  LOBS 

acc  e 

p tab  1 e  for  u 

4) 

MSC  AN 

Cle 

ars  mark 

s  f  r 

om 

atab  ase. 

5) 

SCNDB 

Ac  c 

umu 1 ates 

cou 

n  ts 

of  marke 

6) 

SUX  I 

Set 

a  up  the 

ab  i 

lity  to  extr 

7) 

XSEC 

Use 

s  two  LOBS  t 

o  determine 

us  i 

ng  the  c 

rite 

r  ia 

for  a  va 

8) 

ZEX  I 

Cle 

ars  appr 

op  ia 

te 

nd i c  ator 

fa  i 

Is  valid 

i  ty . 

previous  fix. 


Dx  ,  Dy 
X  e  f  i  x 
Ye  f  i  x 
TempXefix 
T  emoY e  f i x 


REAL; 

REAL; 

REAL, 

REAL; 

REAL; 


<#  x  coordinate  of  estimate  *) 

(#  y  coordinate  of  estimate  *) 

<*  temporary  variable  for  Xefix  *) 
<*  temporary  variable  for  Yefix  *) 


XI.  Y1 

REAL;  ( 

Intx,  Inty 

REAL; 

Bin 

INTEGER; 

B i nNumb  er 

INTEGER; 

Nsame 

INTEGER; 

No  ther 

INTEGER; 

Nmax 

INTEGER; 

Ds  p  x  i 

INTEGER; 

Lobcounter 

INTEGER ; 

Newcoun ter 

INTEGER ; 

Cntfunction 

INTEGER ; 

StaNumb 

INTEGER; 

StationCntr 

INTEGER; 

Node 

BOOLEAN; 

Result 

BOOLEAN; 

Computed 

BOOLEAN; 

An  suier 

BOOLEAN; 

Response 

CHAR; 

DbRecor d 

Lob  type. 

(*  normal  or  degraded  *) 


BEGIN 

WR ITELN ( 

Dx 

Dy 

Xe  f  i  x 
Yefix 
Nsame 
Notber  : 
Nma  x 


Entering  ESTMP ' ) ; 


-  00.  OOC; 
=  00.  OOOi 
=  00.  OOOi 
“  00.  OOOi 
=  0; 


Nmax  : =  Oi  (*  maximum  LOB  counter  ^* ) 

LQADDATA(Lobrec , StaNumb )  i 

FIRSTP ASS( Lobr eci  StaNumb.  Table.  Intio); 

IF  IntPass  >  1  THEN 
BEGIN 

I  NELL  IPSE  ( Xef  i  x .  Yefix.  Answer.  Table). 
TempXefix  :=  Xefix; 

TempYefix  :=  Yefixi 

IF  Answer  =  TRUE  THEN 

BEGIN 

ZEXI(Table,Lobrec  )i 
CLEAR ( Lob rec .  Inti o. Tab  le ) i 

END 

ELSE 

LfiARK (Table.  Lobrec)  i 
WHILE  Index  >  -1  Do 


BEGIN 
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IF  Lobrec  CBin3.  Flag  *  TRUE  THEN 
BEGIN 

Cntf unc  t i on  :  =  Oi 

SCNDB  (Cntf unction,  Index,  Lobrec); 
END 
ELSE 

Lobcounter  : =  Lobcounter  +  1; 

Index  : =  Inde  x  -  1; 

END; 

IF  Nother  =  0  THEN 
BEGIN 

IF  Mode  *  TRUE  THEN 
IF  Nsame  >=  4  THEN 

BEGIN 


Newcaun ter 

:  —  Nsame  -r  Nother; 

IF  Nma  x  >* 

Neucounter  THEN 

BEGIN 

Nma  x  ;  = 

Neucounter; 

Xefix  ; 

=  Intx; 

Yefix  : 

=  Inty; 

Nsame  : 

a  0 ; 

Nother 

:  *  0; 

Neujc  oun 

ter  :  *  0; 

END 

ELSE 

BEGIN 

ZEXI(Table,  Lobrec); 

CLEAR (Lobre c ,  Intlo,  Table); 
Nsame  .  =  0; 

Nother  : =  0; 

Nma  x  :  *  0; 

MSCAN<  Lobre  c ) ; 

END; 


END; 


END 


ELSE 


Newcounter  :=  Nsame  +  Nother* 
IF  Nmax  >=  Neucounter  THEN 


i 


i 

f; 


END; 


BEGIN 


Nma  x 
Xef  i  x 
Yef  i  x 
Nsame 
Noth  er 


=  Newcounter; 
=  I  n  t  x ; 

=  Inty* 

=  0; 

:  =  0; 


END; 


END; 


t 


WR  I TELN  ( 


Leaving  ESTMP  ' )  ; 


I 


END;  <*  PROCEDURE  "ESTMP"  *> 


i 

% 

1 


r.r 

TV 


<  #*******#****#**^*#^******aH*******^***Hfr*********  ****************  *****-*-&■«■** 

<  *************** ****************************  ****************  ****************** 
(  *******  ****  ********  ********  ****************************  ************  ****  ****** 
<*  MAIN  PROGRAM  *) 

BEGIN 

REWRITE(TrailOut) ; 

RESET( Trail  In); 

WR I  TELN  (Tr  a i 1  Out ) ; 

WR ITELN(TrailOut) ; 

WR  ITELN  (Tr  a  i  IQu  t. 

WR  I  TELN  (Tr  a  i  lOut. 

WR  ITELNCTrailOut, 

WR ITELN(TrailOut) 

WR I TELN (TrailQut) 


****  BEGINNING  PROGRAM  TRAILBLAZER 


****  •  ) ; 


Index  : =  0; 

IntPass  1;  ( *  First  pass  through  the  system  *) 


CLEAR ( Lob ,  Int 1 o, Tab le ) ; 


ESTMP(Lob»  IntPass*  Index*  Intlo*  Table); 


yv 


PQ 


WR ITELN (TrailOut); 
WR  ITELN  (TrailOut); 
WR  I  TEL  N(TraiIOut,  ' 
WR  ITELN(Trai  I  Out,  ' 
WR  i  TEL  N(  TrailOut,  ' 
WR  ITELN  (TrailOut); 
WR I TEL  N(TrailQut); 


****  PROGRAM  TRAILBLAZER  COMPLETED  ****'); 


V  END. 


<*  MAIN  PROGRAM  "TRAILBLAZER"  *> 
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*#**  BEGINNING  PROGRAM  TRAILBLAZER  **** 


& 

& 


How  many  stations  will  be  reporting  LOB"s 
on  the  transmitting  object. 

4 

Input  the  Station  that  is  "FIXING". 

1 

Input  the  x-coordinate  of  Station  41 
10.  000 

Input  the  y-coord irate  of  Station  #1 
18.  000 

Input  the  angle;  in  relation  to  true  North/  formed 
by  Station  #1  and  the  possible 
location  (  FIX  )  of  the  transmitting  object 
50.  000 


Input 

th  e 

Station  that 

2 

i  s 

"FIXING 

ll 

Input 

th  a 

x-coord mate 

of 

Station 

42 

20  OOO 

Input 

th  e 

y-coordinate 

of 

Station 

42 

5.  000 


Input  the  angle,  in  relation  to  true  North,  formed 

U  y  ovauiUM  ire:  cjiiu  til  c  pubsiuic 

location  <  FIX  )  of  the  transmitting  object 
33.  000 


Input  the  Station  that  is  "FIXING". 
3 


Input  the  x-coordinate  of  Station  #3 
60.  OOO 

Input  the  y-coordinate  of  Station  #3 
PS  OOO 


Input  the  angle,  in  relation  to  true  North,  formed 

by  Station  #3  and  the  possible 

location  (  FIX  )  of  the  transmitting  object 

-5.  OOO 
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P 

to 


13 


$ 


I 


/V 

VVJ 


& 

W 

& 

fi 


t- 


v.l 


p/ 

.,4  | 


Input  the  Station  that  .is  "FIXING". 


* 


Input  the  x-coordinate  of  Station  #4 
90.  000 


Input  the  y-coordinate  of  Station  #4 
25.  000 

Input  the  angle*  in  relation  to  true  North*  formed 
by  Station  #4  and  the  possible 
location  (  FIX  )  of  the  transmitting  object 
312.  000 


? |  The  two  stations  participating  in  the  intersection 
P  calculations  are; 

Station  #1 
Station  #2 

B 

The  following  information  determines  the  back 
LOB " s  t  o ; 

Station  #3 
Station  #4 


r: 


- Beginning  Of  Inputted  Data - 

— For  Intersecting  Station"*  1  And  2- 


k 


Intersecting  coordinates  for  Station  #1  and  Station  #2  is; 


(50.525*  52.005) 

$ 

i6  X-c  o  ord  i  nate  60.000  of  Station  #3 
^loaded  into  the  system. 

Y-coordinate  25.000  of  Station  #3 
^loaded  into  the  system. 


(V  The  exact  computed  LOB  angle  is  *  -19.334 

C-. 

K 

Difference  between  the  back  LCD  angle  and  ths 
factual  angle  for  Station  #3  -  14.334. 

i 


ty  X-c  o  ord  i  nate  90.000  of  Station  #4 
loaded  into  the  system. 

*  V 

Y-coordinate  25.000  of  Station  #4 
n  loaded  into  the  system. 

D-31 


The  exact  computed  LOB  angle  is 


55.  624 


Difference  between  the  back  LOB  angle  and  the 
actual  angle  for  Station  #4  “  7.624. 


- End  Of  Inputted  Data - - - 

— For  Intersecting  Station"*  1  And  2 - 


Y-: 

& 


The  two  stations  par t ic i pati ng  in  the  intersection 
calculations  are: 

Station  #1 
Station  #3 

The  following  information  determines  the  back 
LOB" s  to: 

Station 

Station  #4 


- . - Beginning  Of  Inputted  Data - — — 

— For  Intersecting  Station"*  1  And  3 — — 


Intersecting  coordinates  for  Station  #1  and  Station  #3  is: 


07.151,  57.564) 


^coordinate  20.000  of  Station  #2 
loaded  into  the  system. 

Y-c o ord i note  5  000  of  Station  #2 
loaded  into  the  system. 

The  exact  computed  LOB  angle  is  *  35,252 

Difference  between  the  back  LOB  angle  and  the 
actual  angle  for  Station  #2  =  2.252. 


/'coordinate  90  000  of  Station  #4 
loaded  into  the  system. 


Y'coordinate  25. GOO  of  Station  #4 


r=w 

SS 


r 


(S 

o 


The  exact  computed  LOB  angle  is  =  -45.249 


Difference  between  the  back  LOB  angle  and  the 
actual  angle  for  Station  #4  *  2.751. 


— End  Of  Inputted  Data — 

- For  Intersecting  Station"*  1  And  3- 


I 


> 


§ 


The  two  stations  par t ic i pat i ng  in  the  intersection 
calculations  are: 

Station  #1 

k-  Station  #4 

k- 

T h  =>  ,  lloujiny  in f oTiTia t i o n  determines  the  back 
LOB " s  to: 

B  Station  #2 

Station  #3 


1 


- Beginning  Of  Inputted  Data - 

-For  Intersecting  Station”*  1  And  - 


r-  N 
»  \ 

£y 


* 

S 


Intersecting  coordinates  for  Station  #1  and  Station  #4  is: 
<55.434,  56.124) 


k:  X-coord  jnate  20.000  of  Station  4)2 
$  loaded  into  the  system. 

-u\  Y-coordinate  5.  000  of  Station  #2 
loaded  into  the  system. 

*  V-’  The  exact  computed  LOB  angle  is  =  34.726 

i  kt* 


»\  Difference  between  the  back  LOB  angle  and  the 
*£\  actual  angle  for  Station  #2  =  1.726. 

ji  X -coordinate  60  000  of  Station  #3 
■r-  loaded  into  the  system. 
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Intersecting  coordinates  for  Station  #2  and  Station  #3  is 


X-coordinate  10.000  of  Station  #1 
loaded  into  the  system. 

Y-cocrdinate  18.000  of  Station  #1 
loaded  into  the  system. 


The  exact  computed  LOB  angle  is  =  46.986 


Difference  between  the  back  LOB  angle  and  the 
actual  angle  for  Station  #1  =  3.014. 


X-coordinate  90.000  of  Station  #4 


loaded  into  the  system. 


Y-coordinate  25.000  of  Station  #4 
loaded  into  the  system. 


i 


>,The  exact  computed  LOD  angle  is  »  -42.  174 


-^Difference  between  the  back  LOB  angle  and  the 
Hactual  angle  for  Station  414  =  5.826. 


Cv  - - - - - - End  Of  Inputted  Da  t  a  - 

JC> — - - For  Intersecting  Station's  2  And  3- 


$ 


& 


The  two 


stat ions 


..  calculat  ions  are  . 

& 


participating 

Station  #2 
Station  44 


in  the  intersection 


The  following 
^;lob"s  to. 


information  determines  the  back 

Station  41 
Station  43 


- Beginning  Of  Inputted  Data - 

— For  Intersecting  Station's  2  And  4 


.Intersecting  coordinates  for  Station  42  and  Station  44  is 

C-} 

<54.024,  57.393) 


^X~coordinate  10.000  of  Station  41 
loaded  into  the  system. 


Y-coordinate  IS. 000  of  Station 
loaded  into  the  system. 


41 


CY 


The  exact  computed  LOB  angle  is 


48.  178 


ft 


Difference  between  the  back 
actual  angle  for  Station  41 


LOB  angle  and  the 
=  1. 822.  d-35 


X-coordinate  60.000  of  Station  #3 
loaded  into  the  system. 

Y-coordinate  29.000  of  Station  #3 
loaded  into  the  system. 


The  exact  computed  LOB  angle  is  ■  -10.452 


Difference  between  the  back  LOB  angle  and  the 
actual  angle  for  Station  #3  ■  5. 452. 


- End  Of  Inputted  Data - 

-For  Intersecting  Station"!  2  And  4 


The  two  stations  participating  in  the  intersection 
calculations  are: 

Station  #3 
Station  #4 

The  following  information  determines  the  back 
LOB" s  to: 

Station  #1 
Station  42 


— — — Beginning  Of  Inputted  Data— - 
-For  Intersecting  Station"!  3  And  4 


Intersecting  coordinates  for  Station  43  and  Station  #4  is 
<57.435.  54.322) 


X-coordinate  10.000  of  Station  41 
loaded  into  the  system. 

Y-coordinate  18.000  of  Station  41 
loaded  into  the  system. 


The  exact  computed  LOB  angle  is  ■  52.558 

D-36 

Difference  between  the  back  LOB  angle  and  the 


«ctu«l  angle  for  Station  #1  ■  2. 998. 


X-coordinate  20.000  of  Station  #2 
load  ad  into  the  system. 

Y-coordinata  9.  000  of  Station  #2 
load  ad  into  tha  system. 


Tha  axact  computed  LOB  angle  is  ■  37. 198 


Difference  between  tha  back  LOB  angla  and  the 
actual  angla  for  Station  #2  *  4. 198. 


> - — - End  Of  Inputted  Data- - - — - - 

- For  Intersecting  Stat i on" s  3  And  4 - 


This  is  the  table  that  shouts  the  relation  of  data 


between 

any  of  the 

stations. 

Intersecting  Stations 

Xintercep t 

Yinter cept 

Supp  ort 

Stat ion 

41 

Station 

41 

0.  000 

0.  000 

0 

Station 

*1 

Station 

42 

90.  525 

52.  005 

2 

Station 

#1 

Station 

43 

57.  151 

57.  564 

0 

Stat ion 

#1 

Station 

44 

55.  434 

56.  124 

2 

Stat ion 

*1 

Station 

45 

0.  000 

0.  000 

0 

Stat ion 

#2 

Station 

41 

50.  525 

52.  005 

2 

Station 

#2 

Station 

42 

0.  000 

0.  000 

0 

Stat ion 

#2 

Station 

43 

56.  793 

61.  656 

0 

.  Station 

#2 

Station 

44 

54.  024 

57.  393 

0 

Stat ion 

*2 

Stat i on 

45 

0.  000 

0.  000 

0 

Station 

#3 

Station 

41 

57.  151 

57.  564 

0 

Stat ion 

#3 

Station 

42 

56.  793 

61.  656 

0 

Station 

#3 

Stati on 

43 

0.  000 

0.  000 

0 

Station 

#3 

Station 

44 

57.  435 

54.  322 

2 

Stat ion 

#3 

Station 

45 

0.  000 

0.  000 

0 

Station 

44 

Stat i on 

41 

55.  434 

56.  124 

2 

Stat ion 

#4 

Stat i on 

42 

54.  024 

57.  393 

0 

Stat ion 

#4 

Station 

43 

57.  435 

54.  322 

0 

Stat ion 

#4 

Station 

44 

0.  000 

0.  000 

0 

*  Station 

44 

Stati on 

45 

0.  000 

0.  000 

0 

Stat ion 

45 

Station 

41 

0.  000 

0.  000 

0 

Stat  i  on 

43 

Stat i on 

42 

0.  000 

0.  000 

0 

Stat  ion 

45 

Station 

43 

0.  000 

0.  000 

0 

Stat i on 

43 

Stati on 

44 

0.  000 

0.  000 

0 

Station 

45 

Stat i on 

45 

0.  000 

0.  000 

0 

_  D-37 

Station  #1  data: 


X-c  oord 

5= 

10.  000 

Y-coord 

as 

18.  000 

Theta 

= 

50.  000 

Alpha 

as 

40.  000 

Reflob 

«c 

TRUE 

Flag 

ss 

FALSE 

Semimajor 

as 

0.  000 

Semiminor 

= 

0.  000 

Distance 

0.  000 

Org  in 

* 

0 

Station  #2  < 

data: 

X-c  oord  = 

20.  000 

Y-coord  = 

5.  000 

Theta  * 

33.  000 

Alpha  * 

57.  OOO 

Reflob  * 

TRUE 

Flag  = 

FALSE 

Semimajor  = 

0.  OOO 

Semi  minor  = 

0.  OOO 

Distance  - 

0.  OOO 

Qrgin  = 

0 

Station  #3  data: 


60.  000 
25.  000 
-5.  000 
95.  000 
TRUE 
TRUE 
O.  OOO 
0.  000 
O.  000 

o 


Station  #4 

d  ata 

X  — r  n  ril'd 

90  OOO 

Y-c  oord 

5= 

25,  OOO 

Theta 

= 

-48.  OOO 

Alpha 

= 

138  OOO 

Reflob 

* 

TRUE 

Flag 

a-. 

FALSE 

Sem l ma  j  o  r 

0.  OOO 

Sem  uninor 

= 

0  OOO 

Distance 

0.  OOO 

Org  i  n 

- 

0 

X-c  o  or  d  - 
Y-coord  = 
Theta  = 

A 1  p  h  d  — 

Reflob 
Flag 

Semimajor  = 
Semiminor  = 
Distance  = 
Qrgin  = 
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Station  #5  data: 

k - 


X-coord 
Y-coord 
tv;  Theta 
lv  Alpha 


a  O.  000 
=  O.  000 
*  0.  000 
a  0.  000 
a  FALSE 
*  FALSE 
O.  000 
O.  000 
O.  000 

o 


Ref  1  ob 

lFla9 

tiSemimajor  =■• 
Semiminor  * 
Distance  = 
{^Orgin  = 


^  From  the  data  submitted  by  each  station  and  the  data 
(vjextracted  from  various  calculations  the  best  "FIX" 
location  of  the  object  attempting  to  transmit  is  at: 


<50.525,  52.005) 


« 

/k‘- _ _ _ _ _ _ _ _ _ 
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APPENDIX  E 

DATA  BASE  ENTRIES  FOR  TRAILBLAZER 


PART  I 

TRAIL-BLAZER 


TABLE  OF  CONTENTS 


1.  Trailblazer  Overview  . 

2-  Utilizes  Analysis  for  HPFIXM  .... 

3.  Attribute  Report  . 

4.  Index  and  Dictionary  Descriptions 


E-3 


TRAILBLAZER 


TRAILBLAZER  Cvervistu 


TRAILBLAZER  PROCESS 

DESCR  IPTIQN: 

TRAILBLAZER  is  a  manned,  ground-based  direction  finding  system. 
The  system  functions  by  intercepting  targets  at  the  master 
control  station  <MCS)  and  providing  at  least  two  DF  bearings 
from  some  combination  of  MCS  and  remote  slave  station  (RSS) 
subsystems.  Successful  deployment  of  the  system  requires 
a  line-of-sight  (LOS)  to  the  target  area  and  to  the  other 
subsystems. 


SYNONYM:  AN/TSG-U4 


Utilizes  Analysis  for  HPFIXM 
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TRAIL-BLAZER 


& 

[;5 

L'1  Utilizes  Structure 


tt 

COUNT 

LEVEL 

NAME 

B 

1 

1.  0  TB_Hpf  i  xm 

r  < 

2 

1.  1  TB_Datset 

b 

3 

1.  2  T  B  _D  i  s  f 

0-* 

4 

1 .  3  TB_Es  tat 

5 

1.  4  TB_Hp  fix 

n 

6 

1.  4.  1 

TB_Aut  c  1 

7 

1.  4.  2 

TB_Dbc  py 

8 

1.  4.  2.  1 

TB_Memr 

YJ 

q 

1.  4.  3 

TE_Estat 

fe 

10 

1.4.  4 

TB  _Es  t  inp 

1  1 

1 .  4.  4  i 

TB_Su  x  i 

12 

1  4.  4  2 

!  B  __ba  x  i 

s 

13 

1.  4.  4.  3 

7B_X  sec 

y\ 

14 

1.  4.  4  3.  i 

7B_I lobs 

15 

1.  4.  4.  3.  1. 

1  TB_Atan2 

© 

16 

1.  4.  4.  3.  2 

TB_Lf  in 

L*t» 

17 

1.  4.  4.  3.  3 

TB_X  c  amp 

IS 

1.  4.  4.  3.  4 

TB_X  prep 

19 

1.  4.  4.  3.  4. 

1  T  B  _C  o  s 

LV 

20 

1.  4.  4  3.  4 

2  TBjSin 

r.-.i 

21 

1.  4.  4  4 

TE_I  ne  1  1 

22 

1.  4.  4.  5 

7E  Lmar  k 

& 

23 

1.  4.  4.  5.  1 

TB_Sc lmk 

B 

24 

1.  4.  4.  6 

TB_Scnd  b 

25 

1.  4.  4.  6.  1 

TB_I pot 

fi¬ 

26 

1.  4.  4.  7 

TB_Zex i 

fe 

27 

1.  4.  4.  B 

TB_Mscan 

20 

1.  4.  5 

TB__Fin  a  1 

29 

1.  4.  5.  1 

TB_Atan2 

1 

30 

1.  4.  5.  2 

TB_C  os 

IP 

31 

1.  4.  5.  3 

TB_F  lobs 

32 

1.  4.  5.  3.  1 

TB_A  tan2 

?■" 

33 

1.4.53.2 

TB_Lf  in 

(V/. 

’v 

34 

1.  4.  5.  4 

TB  _F  pot 

35 

1.  4.  5.  4.  1 

7B_Scnd  b 

36 

1.  4  5  4  1 

1  TB_Ipot 

E; 

37 

1.  4.  5.  4.  2 

TB_Sqrt 

L\ 

3S 

1.  4.  5.  5 

TB _F  sqr  t 

39 

1.  4.  5  6 

TB_Sin 

fv 

f.  ^ 

40 

1.  4.  5  7 

7B_Bcnd  b 

IV  a 

41 

1.  4.  5.  7.  1 

TB_I pot 

42 

1.  4.  5.  3 

TB^Sqrt 

-V 

43 

1.  4.  6 

TB  _Fp  e  a  k 

44 

1.46.  1 

TB  _F lobs 

& 

45 

1.  4.  6.  1.  1 

TB_Atan2 

46 

1.  4.  6.  1.2 

TB_Lf  in 

S' 

v  *■ 

4/ 

1.  4.  6.  2 

TB_Fmneuj 

& 

48 

1.  4.  6.  2.  1 

TB_Fmsam 

/y  ; 

49 

1.  4.  6.  3 

TE_F(nsam 

i  / 

ft 
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TRAILBLAZER 


COUNT  LEVEL  NAME 


50 

1.  4.  6.  4 

TB_Fms  e  t 

51 

1.  4.  6.  5 

TB_Lmar  k 

52 

1.  4.  6.  5.  1 

TB_Sc lmk 

53 

1.  4.  6.  6 

TB  _M  k  a  c  t 

54 

1.  4.  6.  7 

TB  __F  pot 

55 

1.  4.  6.  7.  1 

"”TB_Scndb 

56 

1.  4.  6.  7.  1.  1 

TB_Ip  ot 

57 

1.  4.  6.  7.  2 

TBJSqrt 

58 

1.  4.  6.  S 

TB  Inell 

59 

1.  4.  6.  9 

TB_Zexit 

60 

1.  4.  6.  10. 

TB  _M  s  c  a  n 

61 

1.4.6.  11. 

TB_Zex i 

.  — i 

O  C- 

1.  4.  7 

T'D  _Gr  o  om 

6  J 

i.  4.  a 

TB_Stc  p  y 

a4 

1.  4.  9 

TB  _Mk  usd 

65 

1.  4.  10. 

TB  _Mc  1  r 

66 

1.  4.  11. 

TB_Sve  1 1 

67 

1.  4.  12. 

TB  __Qd  c  p  y 

63 

1.  4.  12.  1 

TB__Mab  i  n 

69 

1.  4.  12.  2 

TB_Mscan 

70 

1.  4.  12.3 

TB  _W  h  b i n 

71 

1.  4.  12.  3.  1 

TB_Lf  in 

72 

1.  4.  12  4 

TB_X  y  2 1  1 

73 

1.5  TE 

_Lamp 

74 

1.6  TB 

__Pc  h  k 

75 

1.  7  TE 

St  ore 

E-6 


TRAIL-BLAZER 


Utilizes  Matrix 


Explanation  of  the  Utilizes  Matrix: 

The  raws  are  input  PROCESS  names;  and  the  columns  are 
PROCESSES  UTILIZED  b  ij  (or  a  SUBPART  of)  the  rows. 


( i ;  j )  value 


mean ing 


Column  j  is  UTILIZED  by  Raw  i 

Column  j  is  a  PART  of  Row  i 

Column  j  is  both  UTILIZED  by;  and  a  PART 


15  TB_Mkusd 
14  TB_Stcpy  -■ 

13  TB_Groom  - 

12  T  0  _F  p  e  a  k - 

11  TB  Final  - 


10  TE_Estmp 
9  TE_Dbcpy  — 

3  TB  _Au  t  c  1 - 

7  TD__Store - 

6  TB  Pc h  k - 


—  / 


S  TB_Lamp  - 

4  TB_Hp fix - 

3  TE_Estat - 

2  TB  _D l s  f  - 

TB  Datset  - 


- / 

/ 


--  / 

/ 


1  TB Hp f i x  m  — — 

2  TB_Hp  fix - 

3  TB_Dbc  p  y - 

4  TB_Estmp  - 

5  TB  Xsec  - 


U 


u  u 
u 


u  u 


u  u  u 


u  u  u 


o 

7 

3 

9 

10 


TB  _  1 1  o  b  s 
T3_Xprep 

TE _ Lma  r k 

TB_Scn  d  b 
TB  Final 


Matrix  empty  for  Rows  11  thru  16  and  Columns  1  thru  15 


E-7 


Li  PJ 


TR AILELAZER 


Utilizes  Matrix 


30  TE_Xcomp 

29  TB_Lf in  - 

28  TB_Ilobs  - 

27  TB  _Ms  can - 

26  TB  Zexi  - 


25  TB _Sc n d b -  / 

24  TB_Lmar  k  - /  ! 

23  TB_Inel  1 - /  !  1 

22  TB_X  sec - /  :  1  ! 

21  TB  Gaxi  -  /  \  \  \  l 


20  TB_J3u  x  i - 

19  TB_Memr  - 

IB  TB_Gdcpy  - / 

17  TB__3ve  1  1 - /  ! 

16  TBJIclr - /  !  ! 


1  TB_Hpfixm - — — — - 

TB_Hp  fix - - - 

TB_Db  c  p  y  - 

TB_Estmp  - 

5  TB  Xsec  - 


U  U  U 


u!  u  u  u  u  u; 


6  TB__I  lobs  - 

7  TB_Xprep  - 

8  TI3_Lmar  k - 

9  TB_Scndb - - - 

10  TB  Final  - 


E-S 


trailelazer 


Utilizes  Matrix 


6  7  8  9  0  1 


3  4  5  a 


11  TB  _F lobs 

12  TB_Fpot  - 

13  TB  _F  peak 

14  TB_Ftnnew 

15  TB  _0  d  c  p  if 

1 a  TB  Wh  bin 


40  T  B  _S  q  r  t - 

39  TB  _Fs  qr t - 

38  TB_F  pot - 

37  T  B  _F  1  o  b  s - 

36  TB_I pot  - 


35  TB_Sclmk 

34  TBJBin  - 

33  TB_Cos  - 

32  TB_Atan2  - 

31  TB_Xpr  e  p - 


1  TB_Hp  f  i  xm 

2  TB_Hp  fix  - 

3  T3_Db  c  p y  - 

4  TB_Estmp 

5  TB  Xs  ec  -■ 


6  T  B  __I  1  a  b  s 

7  TR_Xprep 

8  TBJ_mark 

9  T  B  _S  c  n  d  b 
10  TB  Final 


11  TB_Flobs  - 

12  T  B  _F  p  o  t  - 

13  T3_Fp  eak  - 

14  TB_Fmneui - 

15  T B d  c  p  y  - 


U  U  U 


U  IJ  L 


16  TB  Wh  b  in 


cn  10 


TRAIL-BLAZER 


Ut  i  1  i  z 


es  Matrix 


48  TB_Xy211  - 

47  TB_Whbin  - / 

46  TD _Ma bin  — - - — - — - /  ! 

i  i 
I  I 

45  TB_Ze  x  it - /  !  ! 

44  TB_Mkac  t -  /  !  !  ! 

43  TB_Fmset - /  !  !  !  ! 

42  TB_Fmsam -  /III  i  ! 

41  TB_Finneuj - /  !  !  !  !  '  I 

111*1  It 

1  1  i  i  >  ii 


i 

TB_Hpf  l  mu 

c- 

TB_Hpf l x 

3 

TB_Db  c  p  y 

4 

TB_Estmp 

5 

TB„Xsec  - 

6 

TB_I lob  s 

T 

TB_Xprep 

8 

TB_Lma'-  k 

9 

TB_Scnd  b 

10 

TB„Final 

1  1 

TB_Flobs 

12 

TB  _F  pot  - 

13 

TB_Fpeak 

14 

TB  Fmneui 

15 

TB_0dcpy 

16 

TB_Whb in 

U  U  U  U  U! 
U  ! 


u  u  u: 


-t- 


ts1! 


£ 


TRAILBLA7-ER 
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i 


ft 


ft 


N 


& 


h> 


L 


r< 
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3  Attribute  Report 


REPORT  SPECIFICATION: 

1  tree-level  H= ' TREE  LEVEL '  C0L=14 

2  me thematical-field  H= ' MATHEMAT ICAL  FIELD'  COL=25 


***  No  SYSTEM- PARAMETERS 

TREE  LEVEL  MATHEMATICAL  FIELD 


TB_Atan2 

leaf 

tr  i  gonome try 

T3_Co  s 
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TRAILBLAZER 


1  TB_.  DEQU 

2  TB  SVELL 


1  T B _HP FIX.  MOD  PROCESS 

DESCR  IPTION: 

HPFIX  is  the  main  module  for  processing  a  fix  command 
received  from  either  operator.  It  computes,  displays,  and 
saves  fix  data  as  required  by  the  fix  commands. 


S.  I  w 


'I  XD 


PROCESS 


DESCRIPTION 

P I  XD  is  jumped  co  by  COMRET  to  complete  the  task  of  process;, 
the  FIX  c  ommand. 


G  TB_  ENQU  PROCESS 

DESCR  IPTION: 

Entry  point  for  .ENQU  system  call 

4  TB_.  DEQU  PROCESS 

DESCRIPTION; 

Entry  point  for  .  DEQU  system  call. 


5  TB_.  PTY  PROCESS 

DESCR IPTION: 

Entry  for  .PTY  instruction  call. 


6  TB_.  QUIT  PROCESS 

DESCR  IPTION: 

Entry  point  for  .  QUIT  instruction. 


7  TB_DATSET  PROCESS 

DESCR  IPTION: 

DATSET  is  called  by  F I XD  when  a  fix  has  been  successfully 
computed.  It  generates  a  formatted  table  of  fix  information 
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called  the  DATSET  generated  fix  table. 

8  T  B  _D  I  SF  PROCESS 

DESCR IPTICN: 

DISF  is  called  by  FIXD  to  display  each  successfully  calculated  fix 
and  by  D1SP  when  processing  a  DISPLAY  command,  or  during  display 
change  processing  to  re-display  fix  information. 


?  TB  ESTAT 


PROCESS 


DESCRIPTION- 
T  h  a  function 
messages  for 
messages  are 


of  this  nioilu  le/proc  is  to  queue  error  status 
display  on  the  AN/UYQ-10  and  to  assure  that 
displayed  for  three  seconds,  Any  module  may 


a  i  i 


ESTAT  whenever  display  of  an  error/status  message  is  required 


10  T B _HP FIX  PROCESS 

DESCRIPTION: 

HPFIX  is  called  bu  FIXD  to  compute  a  fix  and  to  provide  all  -  ; 
information  in  a  FIX  DATA  TABLE  containing: 
the  fix  point  latitude 

the  fix  point  longitude 

the  fix  ellipse  orientation  angle  major  axis  relative  to  NORTH 
the  length  in  km  of  the  fix  error  ellipse  semi-major  axis 

the  length  in  km  of  the  fix  error  ellipse  semi-minor  axis 

the  display  unit  number 
the  LOB  display  header 

the  LQB's  involved  in  the  fix  computation. 


11  TB  LAMP 


PROCESS 


DESCR IPTIQN: 

LAMP  is  called  to  program  panel  button  lights.  It  combines  the 
panel  unit  number  and  the  desired  function  into  a  command  word- 
and  sends  the  command  to  the  selected  panel. 


12  TB_PCHK  PROCESS 

DESCR  IPTIQN: 

This  ruuiiny  sun  ply  cheuks  i  u  set?  if  Lilt?  pc*  ye  Lurreiiily  u  ii 
display  has  function  which  depend  on  the  page-  e.  g.  to  see  if 
the  current  page  number  is  one  which  permits  text  editing 


12  TD  STORE 


PROCESS 
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DESCR  IPTION: 

STORE  is  called  by  CQMRET  as  a  result  of  a  STORE  command  and  by 
FIXD  for  each  fix  in  a  multi-fix  request.  STORE  saves  the  late 
set  of  fix  data  in  FIXQ 


14  TB_.  FORK  PROCESS 

DESCRIPTION; 

Entry  point  for  .FORK  system  call. 


1  “  ;  E_AUTCL  PROCESS 

DESCR  IPTION- 

AUTCL  is  called  by  HPFIX  before  computing  the  first  multi-fix. 
and  by  CLEAR  to  process  the  fix  portions  of  the  PURGE  and  DONE 
commands  or  in  response  to  the  ELIM  command  where  no  specific 
fixes  were  specified. 


Id  TBJJBCPY  PROCESS 
DESCR  IPTION: 

DDCPY  is  called  by  HPFIX  to  obtain  a  local  copy  of  that  pari;  of 
the  LOB  database  which  is  displayed  and  not  removed.  The  gioba 
database  referred  to  is  DATAO /DAT A1 . 


17  TB  JESTMP  PROCESS 
DESCR  IPTION: 

ESTflP  is  called  by  HPFIX  to  obtain  a  best  fix  estimate  for 
either  a  single  fix  or  one  of  multiple  fixes  in  either  the 
normal  or  degraded  mode. 


13  TB  FINAL 


PROCESS 


DESCRIPTION: 

FINAL  is  called  by  HPFIX;  after  an  optimized  fix  point  has  been 
determined;  to  calculate  the  parameters  of  the  error  ellipses. 


19  TB_FPEAK  PROCESS 

DESCR IPTION 

FREAK  is  called  by  HPFIX  to  optimize  the  fix  estimate 


E-15 


TRAILBLAZER 


obtained  from  ESTMP. 


20  TE_GRQOM  PROCESS 

DESCRIPTION: 

GROOM  is  called  by  HPFIX  to  initialize  the  local  database  for 
fix  computations  or  to  r e~in i t ia 1 i z e  it  if  the  fix  mode  has  been 
degraded  to  reconsider  rejected  intersections. 


21  TB_MCLR  PROCESS 

DESCR  I  FT  I  ON: 

MCLR  is  a  globally  available  utility  used  to  update  the  header 
text  and  draw  the  page  border  divider  and/or  graphics  for  a  page 
on  the  appropriate  AN/UYQ-10  display  unit. 


22  TB_ODCPY  PROCESS 

DESCR IPTION: 

CDCPY  is  called  by  HPFIX;  after  successful  computation  of  a  fix, 
to  store  all  fix  related  data  in  a  file. 


23  TB_ST CPY  PROCESS 

DESCR IPTION: 

STCPY  is  called  by  HPFIX  to  store  a  local  copy  of  station 
locations  and  screen  geometry  for  use  in  fix  computations 


24  TB_SVELL  PROCESS 

DESCR  IPTION: 

SVELL  is  called  by  HPFIX  to  save  fix  and  related  ellipse 
p  aram  e  ter s . 


25  TB_XY2LL  PROCESS 

DESCR  IPTION: 

XY2LL  is  called  by  HPFIX  following  completion  of  a  successful 
fix  It  takes  X. Y  screen  coordinates  and  performs  a  translation 
and  rotation  to  get  X,Y  coordinates  relative  to  North  being 
directly  vertical  An  inverse  Gnomonic  projection  is  the  done 
to  get  the  latitude  and  longitude  of  the  point. 
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26  TB_MEMR 

DESCR  IPTIGN: 

MEMR  services  memory  requests. 


PROCESS 


I 


SYNONYM:  TB  .  MEMR 
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27  TB_AT AN2  PROCESS 
DESCR  IPTIGN: 

A  TANS  is  <3  globally  available  utility  which  determines  the 
floating  point  radian  angle  whose  tangent  is  X/Y,  X  and  t 
are  double  precision  arguments  supplied  to  the  function 


C- 

IA 


2S  TB_COS  PROCESS 

DESCR  IPTIGN: 

CDS  is  a  globally  available  utility  which.'  passed  an  angle  in 
floating  point  radians,  calculates  its  cosine. 


29  TB_SIN  PROCESS 

DESCR  IPTIQN: 

SIN  is  a  globally  available  utility  which,  given  an  angle  in 
floating  point  radians,  calculates  its  sine. 


£  j 

Wf!  30  7  I  PROCESS 

DESCR  IPTIGN: 

®  SQRT  is  a  globally  available  utility  used  to  approximate  the 

square  root  of  an  input  argument. 
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APPENDIX  F 

USAMS  ALGORITHM  ANALYSIS  SERIES 


Analysis  of  Geographic  Transformation  Algorithms 

JPL  D-181 

DTIC  #ADA  129182 

Dated:  July  9,  1985 

Correlation  Algorithm  Report 
JPL  D-182  UAA-003 

DTIC  0ADA  129181 
Dated:  September  15,  1982 

Applications  of  Correlation  Techniques  for  Battlefield 
Identification  I 

JPL  D-179  UAA-006 

Dated:  June  1984 

Cross-Correlation:  Statistics,  Templating,  and  Doctrine 

JPL  D-184 

DTIC  #ADA  155624 

Dated:  February  29,  1984 

Intelligence  Algorithm  Methodology  I 
JPL  D-183  UAA-004 

DTIC  #ADB  078293 
Dated:  August  15,  1983 

Intelligence  Algorithms  in  Target  Analysis  and  Planning 
(TAP) 

JPL  D-178  UAA-007 

DTIC  #ADB  0924Q2L 

Dated:  November  30,  1984 


7.  Intelligence  Algorithm  Methodology  II:  An 
Intelligence/Electronic  Warfare  (IEW)  Tactical  Sensors  Model 

JPL  D-185  UAA-008 

Dated:  1985 

8.  Intelligence/Electronic  Warfare  (IEW)  Direction-Finding 
and  Fix  Estimation  Analysis  Report 

Volume  1,  Overview 

JPL-180,  Vol.  1  UAA-001 

9.  Intelligence/Electronic  Warfare  (IEW)  Direction-Finding 
and  Fix  Estimation  Analysis  Report 

Volume  3,  GUARDRAIL 

JPL-180,  Vol.  3  UAA-001 

Dated:  December  1985 

10.  A  Non-Standard  Probabilistic  Position-Fixing  Model 

JPL  D-186  UAA-009 

Dated:  June  1985 

11.  A  Collection  of  Area  of  Interest  (AOI)  Algorithms 

JPL  D-171  UAA-011 

Dated:  July  1985 

12.  Power  of  Statistical  Tests  Used  in  Correlation  Techniques 

for  Battlefield  Identification 

JPL  D-2793  UAA-016 

Formerly  Technical  Memorandum  No.  5 

Dated:  August  1985 

13.  Testing  and  Combination  of  Confidence  Ellipses:  A  Geometric 
Analysis 

JPL  D-2782  UAA-013 

Formerly  Technical  Memorandum  No.  2 
Dated;  August  5,  1385 
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14.  Wild  Bearings  Analysis 

JPL  D-2783  UAA-014 

Formerly  Technical  Memorandum  No.  3 
Dated:  July  10,  1985 

15.  Collection  and  Analysis  of  Specific  ELINT  Signal  Parameters 

JPL  D-2781  UAA-012 

Formerly  Technical  Memorandum  No.  1 
Dated:  June  23,  1985 

16.  IEW  Sensor  Error  Budget  for  DF  Fix  Estimations 
Technical  Memorandum  No.  4 

Dated:  August  14,  1985 

17.  Confidence  Ellipse  Research  Software 

JPL  D-2786  UAA-015 

Technical  Memorandum  No.  6 

Dated:  August  8,  1985 

18.  The  Power  of  Statistical  Testa  -  Software 

JPL  D-2788  UAA-017 

Technical  Memorandum  No.  7 
Dated:  December  2,  1985 

19.  Collection  and  Analysis  of  Specific  Elint 
Signal  Parameters:  Final  Report 

JPL  D-2787  UAA-016 

Technical  Memorandum  No.  8 
Dated:  December  9,  1985 

20.  Intelligence/Electronic  Warfare  (IEW)  Direction-Finding 
and  Fix  Estimation  Analysis  Report 

<t  -  i  .  . n  wn  ^  tt  nt  a  orn 

vuuuue  t. )  x iv/i x u jun u t\ 

JPL  D-180,  Vol .  2  UAA-001 

Dated:  December  20,  1985 
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